Formal differential operators, vertex operator algebras 

and zeta- values , II 

Antun Milas 
Department of Mathematics 
University of Arizona, Tucson, AZ 85721 

Abstract 

We introduce certain correlation functions (graded (/-traces) associated to vertex 
operator algebras and superalgebras which we refer to as n-point functions. These 
naturally arise in the studies of representations of Lie algebras of differential oper- 
ators on the circle I22 | -|23 | . |25|. We investigate their properties and consider the 
corresponding graded (/-traces in parallel with the passage from genus to genus 
1 conformal field theory. By using the vertex operator algebra theory we analyze 
in detail correlation functions in some particular cases. We obtain elliptic trans- 
formation properties for (/-traces and the corresponding (/-difference equations. In 
particular, our construction leads to certain correlation functions and (/-difference 
equations investigated by S. Bloch and A. Okounkov 

1 Introduction 

This is a continuation of [25]. In this part we study certain correlation functions built up from 
the iterates of vertex operators introduced in Part I |25|. 

Let V be an arbitrary vertex operator (super)algebra and M a V-module. Sometimes we 
will weakened this property by assuming that M is only a subspace of a IZ-module invariant 
with respect to certain operators. Let Ui € V, i = 1, ...,n, w £ M and w' G M' , where M' is the 
(restricted) dual space of M as defined in |10| . Suppose that Xi's are commuting formal variables 
as in Part I |25| . In the vertex operator algebra theory one usually studies the following formal 
matrix coefficients 

{w' ,Y(u!,xi) ■ ■ -Y(u n ,x n )w), (1.1) 
and the corresponding graded (/-traces (cf. |29|) 

tr| M Jr(ui,a!i) • ■ ■ X{u n ,x n )q m . (1.2) 

Matrix coefficients of the type (II. Ill are studied in 1111 in both the formal and the analytic 
context (meaning that Xi are set to be complex variables). In particular in |12l expressions of 
the form 11. It were used for construction of the genus zero meromorphic conformal field theory. 
On the other hand the graded traces function of the form 11. 21 are related to the genus one 
meromorphic conformal field theory. More precisely, the expression of the form 11.21 give a 
vector in a genus-one conformal block associated to V |29|. 

In Part I |25| we have studied a relationship between various normal orderings and classical Lie 
subalgebras of the Lie algebra of (super)differential operators on the circle. The most interesting 
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representations of T) (unitary representations for instance) can be constructed by using certain 
quadratic operators in terms of free fields |17|. As we already noticed, these quadratic 

operators are closely related to iterates X(Y [u, y]v, x) (cf. |22 | -|21 | . |25|). Therefore it is very 
natural to consider matrix coefficients of the form 

(w',X(Y[ui,yi]vi,xi) ■ ■ ■ X(Y[u n ,y„]v„,x n )w), (1.3) 

which we will call iterated 2n-point function . Also we consider a 2n-point function : 

(u' n+1 , X(ui, xiti)X(vi, xx) ■ ■ ■ X(u n ,x n t n )X(v n ,x n )u n+ i). (1.4) 

The first result in our paper is Theorem 12.11 which relates the matrix coefficients 11.31 and 11.41 . 
In parallel with II. H and 11.21 it is natural to consider the corresponding graded g-traces: 

tr\ M X(Y[a 1 ,y 1 ]v 1 ,x 1 ) ■ ■ ■ X(Y[u n ,y„]v n ,x n )q L( -°\ (1.5) 
Here L(0) = L(0) — Again, (cf. Proposition 12.21 11.51 is closely related to 

tr\iuX(ui,xiti)X(vi,xi) ■ ■ ■ X(u n , x„t n )X(v n , x n )q L( *°\ (1.6) 

Besides the g-traces that we already mentioned, one would like to consider (for reasons that 
will become clear in a moment) an expansion of the q-trace 11.61 in powers of the Xi's: 

tr\ M X(ui,xiti)X(vi,xi) ■ ■ ■ X (u n , x„t n )Y (v„, x n )q L{0) = 

y~](tl\MX(ui,Xitl)X(vi,Xi) ■ ■ ■ X(u n ,X„t n )X(v 
a 

(here we use the multi-index notation x a = a;" 1 • • • x"" ). Especially interesting is the constant 
term with respect to the i-variables, i.e., 

Coeffg.0 , .2.0 ti\MX(ui,xiti)X(vi,xi) ■ ■ ■ X(u„,x„t„)X(v„,x„)q Lm = 

tr| M o(X(ui,a;i^i)X(«i,a:i)) ■ • ■ o(X(u n , x n t n )X(v n , x n ))q L(0) (1.7) 

where o(a) — a(wt(a) — 1) (cf. 29 ). This new formal expression depends on n-parameters: 
ti,...,t n . Also we are interested in 

Coeff^o ^otr| M A(Y"[ui,2/i]«i,:ri) • ■ • X(Y[u 

tr|M0(y[«i,J/l]«l) ■ ■ • o(Y[u n , y n ]v n )q Lm . (1.8) 

Any expression of the form 11.71 or 11.81 (possibly normalized by tr|Mg L ' ') we shall refer to as 
n-point correlation function (or simply, n-point function). This term is widely used in statistical 
physics and random matrix theory where an n-point correlation function depending on yi, ...,y n 
(or ti,...,t n ) can be defined, for example, as an integral (usually normalized by the partition 
function) of a probability density function depending on integration variables sci, x n . Because 
extracting the zero-th term in 11.71 resembles (complex) integration in variables x\,...,x n we 
decided to use this terminology (cf. 3 ). 

In our applications M will be a V-module stable under the Fourier modes of X(tn, xt)X(i)i, x). 
In the case of rational vertex operator algebras with some additional properties (cf. (29p. the 
(/-trace 11.21 rise to a doubly periodic functions (in this case even more is true; the vector space 
of all characters of V-modules is invariant under SL(2, Z) action). However, to prove ellipticity 
one needs much weaker assumptions. 

Here is a short overview: 

1 2n refers to the number of formal (or complex) variables 
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(i) In Section 2 we derive a precise relationship between 11.31 1 and 11.41 . We also relate the 
corresponding g-traces 11.51 and 11.61 and the corresponding n-point functions 11.71 and 

(ii) In Section 3 we apply our results from Section 2 in the case of the infinite-wedge vertex 
operator (super)algebra T and its charge m subspace T m - It is well-known that these 
subspaces are D-modules. We compute q-traces 11.51 and 11.61 in the most interesting 
case (when the vectors m and w, are chosen to be the fermionic generators) and the 
corresponding n-point functions. In particular, the n-point correlation functions 11.81 rise 
to the n-point functions introduced by Bloch and Okounkov [3], |27|. The rest of Section 3 
is devoted to studies of the g-difference equations satisfied by these n-point functions. We 
obtain several explicit formulas for 1- and 2-point functions by using different methods. 

(iii) In Section 4, along the lines of Section 3, we consider a fermionic vertex operator super- 
algebra and the corresponding n-point functions associated with it. This case rises to a 
representation of an orthogonal Lie algebra of differential operators on the circle T>~ . 

(iv) In Section 5, we consider a vertex operator algebra associated to a free boson and the cor- 
responding n-point functions associated to the bosonic generators. This case was studied 
in |25| in connection with a symplectic Lie algebra of differential operators on the circle 
T> + [3], |25| and the zeta-regularization procedure. 

(v) In Appendix A we prove certain elementary lemmas necessary to deal with q-difference 
equations in sections 3,4 and 5. Finally, in Appendix B we gave a different proof of the 
so-called "recursion formula" for the q-graded traces originally due to Zhu 29 . 

n.b. This paper was originally written in 2000 and it is a part of authors Ph.D. thesis. In the 
meantime we have noticed that Okounkov |27| . Miyamoto |26l and Dong, Mason and Nagatomo 
7 obtained some interesting results related to our paper. 

2 Correlation functions 

As in Part I |25| we will use formal calculus as developed in |11| . We denote by x, y, t, Xi, U and 
yi, etc., commuting formal variables and we take the liberty of using the same notation when 
formal variables are replaced by complex numbers. From the context it should be clear whether 
the variables are formal or complex. 

2.1 Normal ordering procedure 

Let ( V, Y, 1, w) be a vertex operator algebra ( see |11| for the definition ) and u,v G V. From the 
Jacobi identity for vertex operator algebras it follows that 



Y{u{ 



l)v,x) = \Y(u,x)Y{v,x)\ := Y~(u, x)Y(v, x) + Y(v,x)Y + (u,x), 



(2.1) 



where Y + (v,x) — v(n)x~ n ~ 

n>0 

part of Y(v,x). Also, we define 



is the singular and Y (v,x) — v(n) 

n<0 



is the regular 



>x 



;Y(u,x)Y(v,y); = Y (u,x)Y(v,y) + Y(v,y)Y + (u,x). 
Let us recall the definition of the A-operator as in the first part |25| : 

X(u,x) = Y(x L(0) u,x), 
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where u is a homogeneous vector. We extend this definition by the linearity for every u € V. 
The X-operator also admits a splitting into the regular and singular part, but then X~ (u, x) 7^ 
Y~ (x L ^u, x). Thus we have two different normal orderings for X operators 

l'X(u,x)X(v,x)l' = X~(u,x)X{v,x) +X(v,x)X + (u,x) 

and 

\X(u,x)X(v,x); = Y~(x L[0) u,x)X(v,x) + X(v,x)Y + (x L(0) u,x). 
It is more convenient to work with • J instead of *' *' . If we suppose that 

u,v€Vi and [v(0), u(n)] = 0, 

for every n, then 

l'X(u,x)X(v,x);' = I ^ii(n)a;" n \ X(v, x) + X(v, x) I J2 u ( n ) x 

[n<0 I Ln>0 

^ u(n)aT n I X(v, x) + u(Q)X(v, x) + X(v, x) I ^ 

n<0 J ln>0 



= < x u{n)x n 1 > k) + A(«, x):r < u(r 

I ra<0 J [n>0 



= xY~(u,x)X(v,x) + X(v,x)xY + (u,x) = ;X(«,a;)X(w,x):, (2.2) 

meaning that there is no ambiguity which normal ordering we are using. This case arises when 
we study bosons and fermions 2 . 

Proposition 2.1 Suppose that 

u(n)v = c„,„<5 niWt(tl ) +wt („)_il, (2.3) 

where c«,„ 6 C. 27ien 



fa) 

and 

(b) 



lX(u,xi)X(v,x 2 ); = ;X(u,s 2 )X(«,a;i): 



:'jf(tt,a5i)X(t;,a;2):' = :X(u, £i)X(«, x 2 ) ; . 
We will not use this result in the rest of the paper so we leave the proof to the reader. 

2.2 Iterated 2n— point functions 

In what follows we shall always use the following binomial expansion conventions. An expression 
[x- y ) k ' wnere x an( A V are formal variables, is understood to be expanded (by the binomial 
theorem) in non-negative powers of y. Note that the order of variables is important. Also we 
allow x or y (but not both) to be a complex number. 

2 All results about normal ordering hold with minor modifications for the vertex operator superal- 
gebras. In particular, the formula 12.11 has to be replaced by • Y (u, x)Y (v, x) • := Y~ (u, x)Y (v, x) + 

(-1)pWpWy(v,x)Y+(u,x). 
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As in 10 we will be considering formal series of the form 

g(xi,...,x n ) 



£C{[xf L ,...,x±% (2.4) 



Yl x 7 Y[( x j ~ x kT 

i= 1 j<k 



where g(x 1 , x n ) € C[x 1 , x n ]. 

It is well-known (cf. HOI, EE]) that 

(u' n+1 ,Y(u-i,xi) ■ ■ ■ Y(u ni x n )u n+1 ) = f(xi,..., x n ), 

for some f(xi,...,x n ) of the form 12.41 1. After we replace the formal variables with complex 
variables 12. II converges to a rational function inside the domain \xi\ > • • • > \x n \ > 0. 

Let pi, Sj,pi,j, rtj, Si t j and be natural numbers and g(xi,tj)™j= 1 £ C[xi, tj]", =1 . Then 

f(xi,t s )? J=1 = (2.5) 

_ 9( X i> tj)™j = l 

n n 

j^J {ti 1) J^J X j J^J (tiXi tj Xj ) ''^ (t {Xi Xj ) {Xi tjXj^j ' J (xj Xj ) ,:? 

is a well-defined element of C^a:^ 1 , x^ 1 , t^ 1 , where we expand the denominator by 
using binomial expansion. Viewed as a rational function 12.51 1 has poles on the divisors: Xi = Xj, 
Xi = 0, U = 0, ti = 1, Xi = - J ^ 2 - , Xi = tjXj and aij = Xjt~ Y . It is also clear that 

{u'„ +1 ,X(ui,xiti)X(vi,xi) ■ ■ ■ X(u n ,x n t n )X(v n ,x n )u n+1 ) = g(xi,tj)i-i, (2.6) 

where f(xi, x n , ti, t n ) is of the form 12.51 . We may assume = ... = ti n for some 
1 < i\ < . . . < ik < n; in this case 12.51 is still well-defined. 

It is important to notice that in 12.51 we cannot perform the substitution ti — e Vi , i = 1, n, 
since the binomial expansion convention is not applicable to the expansion r^Zrp' • R* a th er we 
use a different convention. As before ^ e y ix ^ e vj x .^k are expanded by using binomial expansion, 
and factors of the form 

1 



( e *_l)fc' 

k g N stand for the formal multiplicative inverse of (e x — l) k , i.e. 

ec((x)), 



(l + (f + £ + ...))* 



where we use binomial expansion with respect to !? + %" + •••■ Sometimes we will view 7^7—^ 



as a meromorphic function in y. With these conventions we show 
Proposition 2.2 /(a;^, e Hj )™ =1 J=1 is a well-defined element of 

C(( JW ,...,y„))[[xf\...,a£ 1 ]]. 
Proof: It is enough to show that in the expansion 
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(we use a multi index notation a = (ai, a„)) for each a, f a (e y ')™= 1 is a well-defined element 
of C((2/i, ...,£/„)). 

CoefiV yv ;t ' J_ 



]^(e M, a;i - e Vj x-j) Pi - j (e v< iSi - Xj) Ti - j (xi - e yj Xj) Si - j (xi - Xjf 
Coeff^i (CoefI> 2 (...Coefl>„ 



ff(^,e^')r, J= i 



Joeff 



-)) 



= Coeff^! (CoefF^ (...Coeff^-, 

-)), (2-7) 



]^[(e yi a;i — e Vj Xj) Pi - j (e Vi Xi — Xj) ri - j (xi — e Vj Xj) aiJ (xi — x : 

i<3 

where h is some Laurent polynomial. By induction it follows that (12.711 is an element of 
C[[ yi ,...,y n ]]. ThuB/„(e«")£=iGC(( Vl ,...,y„)). I 

Definition 2.1 Let m, Vj G V, for i = 1, n + 1, j = 1, n and tt^ +1 G V . We define an 
iterated 2n-point function as a formal matrix coefficient: 

(u' n+1 ,X(Y[u 1 ,y 1 ]v 1 ,xi) ■ ■ ■ X(Y[u n ,y n ]v n ,x„)u n+1 ), (2.8) 

where 

Y[u, y] = e wtiu)y Y(u, e v - 1). (2.9) 

Vertex operators of the form 12.91 1 have been introduced by Y. Zhu in |29l . A quadruple 
(V, Y[ , ], l,a>), where ui — ui — ^j, is a vertex operator algebra isomorphic to (V, Y{ , ), 1, w). 
Operators of the type X(Y[u, y]v, x) were introduced in |22|. |23| in connection with the Riemann 
(-function (cf. 

2.3 Associativity for iterated 2n— point functions 

The aim is to relate 12.81 and 12.61 . Let us consider the simplest case when n = 1. We know (cf. 
) that 

X(Y[u,2/]u,a;) = ;X(u,e v x)X(v,x)l + X(Y + [u,y]v,x), 

where 

Y + [u,y] = Y.< n ^ eV 1 

n>0 

and the normal ordering • • for X -operators is the one introduced in Section f2. II Hence 
X(Y[u,y\v,x x ) = :X{u,e y x)X(v,x): +X(Y+[u,y]v,x) 

= ; X(u, e>x)X{v, x) I + ^(eT t(u) S^+j ■ (2-10) 

On the other hand 

X(u,tx)X(v,x) = 

= lX(u,tx)X(v,x)l + [Y + ((tx) L(0) u,tx),X(v,x)] (2.11) 



6 



Now we work out the second term in 12. 1H . 
[Y + {{tx) L(0) u,tx),X(v,x)] = 

= (tx) v * (u) x vn(v) lim Sing x , \Y + (u,tx x ),Y{v,x)] = 

x\ — >x 

= (tx) vrt ^x wt ^ lim Sing Res X(j x~ 1 5 ( — — — ) Y (Y (u, xo)v, x) 

Xl-tX 1 \ X J 

= (tx) wt(u) x wt(l,) lim Sing^ Y x- 1 e- x °^^s( — ')Y(u({)v,x) 

i>0 v y 

= jwtw^wtw+wtM ^ g ' y («(«>,») = 
= ^ t wt(«) a , w t(«)+wtw-i-i — L_y ( u (i) V) x) = 



i>0 

iwt(u) 

i>0 ^ ' 



where Sing stands for the regular part with respect to x\ and lim stands for the formal 

x\ — >x 

substitution xim (cf. 
It is obvious also that 

(«2, J-^Cw, tx)X(v, x) IU2) G C[2, a;,t — ,x — ], 

(u' 2 , ;X(u,e v x)X(v,x)lu 2 ) G C[e v ,x,e~ v ,x~ 1 ], (2.13) 

for u 2 G M and u' 2 G M'. 

If we combine l|2.1()|l . I|2.11^ and 1)2.120 we obtain 

Proposition 2.3 There exist a polynomial g such that 

{u' 2 ,X(u 1 ,tx)X(v 1 ,x)u 2 ) = pjjjr^p 

/or some p,r,s G N, <7(x, 3/) G C[at, j/]. Moreover, 

(uLX(Y[ui,y]vi,x)u2) = 9{e v x,x) 
\ 2, 7 (e») r as s (e* - 1)p 

Jn particular, when x and y are specialized to be complex variables 

{u' 2 ,X(Y[ui,y]vi,x)u2} 

converges to a holomorphic function inside \x\ > 0, < \y\ < 27r. 

Proof: We may assume that all vectors are homogeneous. Relations flTEfl . i2TH and <2"T2l 
imply that 

(U 2 , X(Y \Ul,y \vi,X)U2) = "7 ; ; -r— . 

,yj ' ; ' (ef) r x s (ef - 1)p 
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The second statement follows from the above discussion and the fact that eV 1 _ 1 , as a complex 
function, has a Laurent expansion i + •• • which converges inside < \y\ < 2tv. This expansion 
coincides with the expansion 

1 1 1 e y - 1 - v 

y +... 



e V-l y(1 + £i^i) y y 

exhibited according to our conventions. 

I 

Let us consider the general case (n > 2). For every n G N we define domains 

Qn.n — \0^1 ? ••") Xn, , ^n) £ C 

|*l£Dl| > \xi\ > \t 2 X 2 \ > \x 2 \ > .... > \t n -iX n -l\ > \x n \ > 0}. (2-14) 

It is easy to see that 12.51 is analytic inside f2 n ,n- Note that |tj| > 1, holds for every i = 1, n and 
(x, t) := [x\, x n ,ti, t n ) £ il n ,n- We have the following theorem (that we call associativity) 

Theorem 2.1 There exist a formal series f(xi,tj)™j= 1 of the form i2.5t such that 

(u' n+ i, X(ui,X!ti)X (vi,xi) ■ ■ ■ X(u n ,x n t n )X(v n ,x„)u n+1 ) = f(x i ,tj)" j ^ 1 . (2-15) 

and 

(u' n+1 X(Y[ui,yi]vi,Xi) ■ ■ ■ X(Y[u n ,y n ]v n ,x n )u„+i) = f(xi,e Vl )" j=1 , (2.16) 

where f(xi,e Vj )"j =1 is a series viewed according to our expansion conventions. In particular, as 
complex functions, 12.15\) and 12. 1 61) have the same meromorphic continuation onto C 2n if we 
replace e Vi by ti. 

Proof: We know that 

(u' n+1 ,X(ui,x 1 t 1 )X(v 1 ,xi) ■ ■ ■ X(u„,x n t n )X(v n ,x n )u n+1 ) = (2-17) 

{u'n + l, {'.X(lll,Xltl)X(vi,Xl) I + Fl(tl,Xl)) ■ ■ ■ 

■ ■ ■ {\X(u\,x\t\)X(vY,x\) I + F n (t n ,x„))u n+ i), 

where 

Fi(ti, x„) = [Y + {{Ux,) L{0) u,tx),X{v,x,)] 

On the other hand 

(u' n+1 ,X(Y[u 1 ,y 1 ]v 1 ,x 1 ) ■ ■ ■ X{Y[u„,y n ]v„,x n )u n+ i) = (2-18) 
= {u' n+1 ,( : X(u 1 ,e^x 1 )X(v 1 ,x 1 ) : +f 1 (e y \x 1 ))... 
. . . (lX (u n , e Vn x n )X (v n , x n ) I + f n (e Vn ,x„))u n+ i), 

where 

fi{e y ',Xi) = X(Y + [ui,yi}vi,Xi), 

and the (operator valued) series on the right hand side are obtained according to our expansion 
convention. 

Because 12.151 1 is a sum of several terms of the form 

I .wtO) w w x .\ 

(«n+i, ;X(ui,xiti)X{yi,xi); 3 — ; ' 31 ■■■ ;X(u n ,x n t n )X(v n ,x n )\u n+1 ), (2.19) 

Kyi L ) 
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for some m > 0, the only thing we have to prove is that 

(u' n+1 , ;X(ui, xih)X(vi,xi) I ■■■ X(w,Xj) ■ ■ ■ iX(u n ,x n t„)X(v„,x n ) \u n+1 ), (2.20) 

is an expansion of a meromorphic function with (possible) poles at Xi = 0, t% = 0, Xi = Xj, 
tiXi = Xj, i ^ j and tiXi = tjXj. If this is the case then I2.20H does not contribute with 
additional poles at t j = 1 so we are allowed to make a substitution ti — e Vi in 12.2011 and use the 
same expansion conventions (which are applicable !) and obtain 

{u'n+it ;X(u 1 ,x 1 e vl )X(v 1 ,x 1 ) I ■■■ X(w,Xj) ■ ■ ■ \X(u n ,x n e Vn )X(v n , x„) ;u„ + i). (2-21) 

But we know that 12. 181 is a sum of several terms like 

i e vi vt ( w ) X(w Xi) 
(u n +i, ;X{u 1 ,x 1 e yi )X(v 1 ,x 1 )l - ^ lX{u n ,x n e Vn )X(v n ,x n )lu n+1 ), 

(2.22) 

so we have the proof. 

Claim: The expression 12.191 is an expansion of a rational function with no poles at ti = l. 
Proof: If we use 

\X(ui,Xiti)X(vi,Xi)l = X~(ui, XiU)X{vi,Xi) + X(vi,Xi)X + (ui,XiU), 
for i = 1 we get 

("n+l, '.X(U1, Xltl)X(vi, Xl) I ■ ■ ■ X(w,Xj) ■ ■ ■ lX(u n ,X„t„)X(v n , X„)lU n + l) = 

= (u' n+1 , X + (in, xiti)X (vi, xi) I ■ ■ ■ X(w,Xj) ■ ■ ■ \ X(u n ,x n t n )X(v n ,x n )\u n+ i) + 
{u' n+1 , X(vi,xi)X~ (ui,xiti) ■ ■ ■ lX(u n ,x n t n )X(v n ,x n )lu n+ i). (2.23) 

After we apply the same formula, for every i, we obtain several summons. A generic term is 

{u' n+1 ,X(v 1 ,x 1 )X' (ui^^!) ■ ■ ■ X(wj,Xj) ■ ■ ■ X + (u„,x„t„)X(v n , x„)u„ +1 ). (2.24) 

It is important to notice that X + (ui, tiXi) is always to the right of X(vi,Xi) and X~ (ui,t%Xi) is 
to the left of X(vi, Xi). 

Now, by using the commutator formula we can move X + (ui,tiXi) all the way to the right 
(such that it acts on w„+i) and X~ (ui,Xi) all the way to the left. This procedure will produce 
poles only at wanted places, i.e., no new poles at ti = 1. After we apply X + (ui,tiXi) to w n +i we 
obtain only finitely many terms. The same thing happens with X~ (ui,Xi). We can repeat the 
same procedure for remaining terms. This is a finite procedure and at the end we end up with a 
rational function with no poles at ti = 1. 

The proof follows. 

I 

We define an n-point function to be a formal matrix coefficient: 

(u' n+1 ,o(Y[ui,yi]vi) ■ ■ ■ o(Y[u n ,y„]v„)u„ + i), (2.25) 

where o(a) — a(wt(a) — 1), i.e., we extract the zeroth coefficient in l|2.8|l with respect to the 
expansion in terms of Xi's. In parallel with the previous construction we will consider: 

(u' n+1 ,o(X(u 1 ,t 1 xi)X(vi,x 1 )) ■ ■ ■ o(X(u n ,t n x n )X(v n ,x n ))u n+ i), (2.26) 
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where 

o(X(u, tx)X (v, x)) := Kes x x~ 1 X(u,tx)X(v,x). 

For n > 1, let 

n n , = {(ti,...,*n) G C n : > l,i= l,...,n}. 

Then we have 

Corollary 2.1 T/ie expression 

/(«!,..-,*»)= (2-27) 
= (uJj + i,o(X(ui,tia;i)X(«i,a:i)) • ■ ■ o(X(ti n , t n x„)JT («„, a; n ))?i n+ i}. 

converges to a holomorphic function inside Q n ,o ond it has a meromorphic extension to C™ with 
poles at ti = 1 and ti = 0. 

g(yi,---,y n ) = ("n+i,o(F[ui,yi]«i)---o(F[u„,j/ n ]t; n )?i n+ i) (2.28) 

converges inside < < 27r. Moreover, as meromorphic functions, 

g(yi, -,Vn) = f(h, ...,i n )|t 1 =ew,...,t„= e v«, (2.29) 

Proof: From the Proposition 12. 21 we know that 12.281 is well-defined. Here we make the following 
additional observation: 



J^J(tiXi Xj) *'l (jVi tjXj) %, i{Xi Xj) 
i<j 



is a Laurent polynomial in ti, i = 1, n. Because of 12. (il . it follows that l|2.27^ converges inside 
t;| > 1. Similarly, 



Coeffj.o...;,.!) 



J^J(c Xi •Ej') ^ ^Xi e x j*) (x$ 

i<3 



is (the same) Laurent polynomial in ti — e Vi , i = 1, ...n and it converges inside < \yi\ < 2n. 
Now, Theorem 12.11 implies the last statement in the corollary. 



Remark 2.1 In Proposition 12. 31 we did not use the associativity (as stated in (!()]). One might 
think that there is an easier proof (cf. [22 ] -|2T ] 1: 

X(u,tx)X(v,x) = (tx) wUu) x wUv) Y(u,tx)Y(v,x) ~ 
(hx^^x^YiY^tx - x)v,x) ~ 
{t lX y t{u) X(x- L{<>) Y(u,tx - x)x L{0) v,x) ~ 

X{t H0) Y{u,t- l)v,x). (2.30) 

Even though 12.31)11 implies the "right" answer (after we substitute e y for t), the proof is not 
rigorous ! Note that the term X(t L ^°'Y (u, t — l)v, x) is a non-rigorous expression (this can be 
seen if we extract zeroth coefficient of i) compared with the left hand side which is well-defined. 
However, if we use matrix coefficients this construction can be made completely rigorous. 
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2.4 q— traces for the iterated 2n— point functions 

In the previous section we showed that two correlation functions are related by a simple substi- 
tution f Theorem 12. H . Here we extend the same result for the corresponding q-traces. Let us 
consider (cf. |29|) a formal series (g-trace): 

Xj\uX{ai,Xi) ■ ■ ■ X(a n ,x n )q L(0) € q h C[[xf\ xt\q}], (2.31) 

tti 6 V, i = 1, ...,n, where M is a V -module which admits N-grading, i.e., there is ft € C such 
that Spec(L(0)) £ h + Z>o- For m, Vi € V, i = 1, n we define a g-trace for an iterated 2n-point 
function as a formal expression 

tr|jwrX(y[t»i,yi]wi,a:i) • ■ • X{Y[u n ,y n ]v n ,x n )q L(0) . (2.32) 

Later q = e 2,r ' T , where t € HI (H is the upper half- plane). 

It is a novelty, comparing to 12.3111 , to consider correlation functions of the form 12.321 . Now 
the vector at is replaced by the formal expression Y[m,y]vi. Thus one can think of 12.321 as a 
g-trace attached to a certain completion of V. 

We will work intensively with another related g-trace: 

tr\ M X(u 1 ,t 1 x 1 )X(v 1 , xi) ■ ■ ■ X(u n ,t n x n )X(u n ,x n )q L<0) . (2.33) 
By the definition: 

tr\ M X{u 1 ,tixi)X(v 1 ,Xi)---X{u n ,t n x n )X(v n ,x n )q H0) = (2.34) 

^ {w'j, X(ui,tiXi)X(vi , Xi) ■ ■ ■ X(u n ,t n x n )X(v n ,x n )wi)q wt( - Wi \ 
iei 

where {ujijie/ is some (or any) homogenuous basis for M such that {w' iy Wi) — 1. We know that 

{w'i, X(ui,tiXi)X(vi,xi) ■ ■ ■ X(u n ,t n x„)X(v n ,x n )wi} 
converges to an analytic function in the domain 

|tl!Cl| > > ... > \t n X„\ > \x n \ > 0, 

but a priori we do not know whether (12.331 has a domain of convergence. To prove the conver- 
gence it suffices (according to |29| 'l to check that for every v 6 V the formal 1-point function 

tl\ M X(v,x)q Lm 

converges to a holomorphic function inside |g| < 1. This statement is very strong and it depends 
on the internal structure of the vertex operator algebra (C2-condition for instance). If this is the 
case then one can show that 12.341 converges to a holomorphic function inside 

1 > |tia!i| > > . . . > \t„x n \ > \x n \ > \q\, (2.35) 

and it has a meromorphic extension to (C x ) 2n . 

In our approach 1-point function does not play a prominent role (actually we start with 
2-point functions). Also we do not assume rationality, C2-condition etc.. Therefore we need 
to develop a right recursion procedure for expressing g-traces of iterated 2n — point functions in 
terms of g-traces of certain (2n — 2)-point functions. Then, by using this result we prove the 
convergence of the iterated 2n-point. We shall adopt this approach at the very end. 

From now on we will assume that 12.31 holds for all pairs of vectors Ui,Vj 6 V, i = 1, ...,n. 
The following fact follows immediately from Theorem 12. II 
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Corollary 2.2 Let us denote 12. 33d by f(t,x,q) and 12.32(1 by g(y,x,q). Then 

g(y,x,q) = f(e y ,x,q). (2.36) 

If we assume that x, y and t are complex variables and q = e 2 ^' T then 12. 360 holds, provided that 
both sides are absolutely convergent in a certain domain. 

The previous corollary is not so useful since we do not have any information about poles. 
In the case n — 1 we have the following description |29| (for the notation and a different proof 
see Appendix A) . 

Theorem 2.2 Formally, 

tr\ M X(Y[u,y]v,x) = (2.37) 
^ p m +i(y,q)tr\ M X(u[m\v,x)q L{ - 0) + tT\ M o(u)o(v)q Lm , 

m>l 

where p m +i(y) (defined in the Appendix B) are considered as elements from C((y))[[g]]. 

Remark 2.2 All results obtained in the previous section carry out (in a straightforward manner) 
for vertex operator superalgebras (cf. |16p. Because of the ^Z-grading, correlation functions are 
multivalued functions and therefore the statements in Theorem 12.11 and Corollary 12.21 are am- 
biguous. If we assume that p(v,i) = pivi) (where e is the sign), in 12.33H possible multivaluedness 
stems only from the term 

y/ti 1 ...ti k , 

for some 1 < i, < ... < ik < n. 

3 q— traces and T>\ charged fermions 
3.1 Charged free fermions 

In this section we exploit our general consideration in the case of the vertex operator superalgebra 
stemming from the pair of charged fermions f [14]. [13p. Also, we discuss an equivalent bosonic 
construction. 

First we introduce free (charged) fermions. We consider a Lie superalgebra generated by %j) n 
and ipmi n,m G Z + 1/2, and 1 with commutation relations 

bPn, 1pm] = <Wn,0- [C 1pm] = bPn, ^m] = 0. 

We denote the corresponding Fock space by T ( |13|') spanned by elements 

if>*(-n k ) ■ ■ ■ tp*(-ni)ip(-mi) ■ ■ ■ ip(-mi)l, (3.1) 

where rii, mj € N + \ and rife > . . . > n\ > 1/2, mi > . . . > mi > i. 
Let us define the following vertex operators (fermionic fields): 

nv>(- 1/2)1, x) = y, Tp( n + m*-"- 1 

nez 

and 

Y(ip* (-1/2)1, x) = y, i>*^ + im x - n -\ 
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We write shorthand ip — 1/2)1 and if)* = ip*(— 1/2)1. Since the fermionic fields are local 
and generate the space T we can equip T with a structure of vertex operator superalgebra (see 
HH) with the conformal vector 

w = ^V*(-3/2)V(- 1/2)1 + ^(-3/2)V*(- 1/2)1. 
This vertex operator superalgebra has the charge decomposition 

•F — -F' n , 

TlgZ 

with respect to the operator 

o{;X{i>,x)X{i>* ,x)l). 

Then JT is a vertex operator algebra and every T n is an .T-o-module. 

The projective representation of the Lie algebra T> of differential operators on C x (see |18|. 
, studied in can be interpreted in the language of the vertex operator algebras in the 
following way. 

Consider the following vectors: 

^(-i-l/2)V>*(-i- 1/2)1, (3.2) 

i, j G N. Then T> has a projective representation in terms of the Fourier coefficients of the vectors 
f3~2l (see [5J). Let V be the Cartan subalgebra of V (see HHJ) spanned by L^ r) = (t^j)"", r> 0. 

We define an action of T>o in terms of generating functions (we choose a particular lifting) in 
the following way: 

where L^ r '(0) is the coefficient of x°y r in 

: X^,e v x)X(r,x)'.+-^j. (3.3) 

According to [5], I22| and [55 (13.311 corresponds to a new normal ordering. The operator 
(0) introduced in coincide with D r for every r G N. It is easy to see that operators L^' (0) 
act semisimply on T and therefore we can define a generalized character of T (cf. |0]): 

tri^n«f W(o) - < 3 - 4 ) 

i=0 

where (ft = e 2,rlTi , for i > 1 where ri G H, and qo, Ti, i > 2 are formal variables. Notice that 
L' ^(0) is not associated to any particular element of T>. 
It is important to notice that 13.31 1 can be written as 

X(Yty,y]r,x). (3.5) 

Remark 3.1 One can derive a closed formula for the commutator 

[X(Y[1>, Vl ]i>, xi),X(Y[i>, V 2W,x 2 )], (3.6) 

by using the Jacobi identity from |25[ or 13.31 . We leave this (non-trivial) exercise to the reader. 
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Then 



trkn«r (o) =n^-T +i ' i/2) 

i=0 i=l 

n a + qoq r i/2 i { r i/2)2 ■ ■ oa + ^r +i/ v n+1/2)2 • ■ o- (3.7) 

n>0 

This corresponds to the formula (1.12) in where the same generating function is denoted 
by n(ro, ti, ...). The infinite product has an expansion as a Laurent series of go- We denote 
by V(t\,T2-..) the coefficient of gf} of the generating function 13.71 1. Q. and V are quasi-modular 
forms (for the definition of quasi-modularity see 0, |25|~) of the weight and — | respectively. 

3.2 Boson-fermion correspondence 

Let I be a rank one lattice with a generator a such that < a, a >= 1. We denote by L the 
non-trivial central extension of L 

1 -> {±1} -> £ L -> 0, 

such that a& = (-l) <s,!,> 6a, a,b e L. 
As in [fj] (or |13p we equip the space 

Af(l) ® C[£], (3.8) 
with a structure of vertex operator superalgebra. Then 



ex P ( E ^^s" j exp ( ^ 

\n>0 / \n>0 



X^ 1 , x) = exp ( V ^^x" ] exp ( V a^z^, (3.9) 



where a = a. The boson-fermion correspondence (cf. |13|) is a vertex operator superalgebra 
isomorphism 

* : T V L , (3.10) 

with the following properties: 

ip(-n + 1/2) ■ ■ ■ ^(-1/2)1 i — > e na 
iP*(-n + 1/2) • ■ ■ i>*{- 1/2)1 h-> e"" a , 

c, _> i a (-l) a l. (3.11) 

In particular, ^(-1/2)1 h-» e Q and V>*(-1/2)1 h-> e" Q . 
Now 113.511 can be expressed in the following form: 

X(Y[ip,y]r,x) = 9- 1 X(y[e a ,y]e- e ',x)9, (3.12) 

since 

y(Y{^,y]r) = Y{e a ,y}e- a , 

which is a consequence of Proposition 4.1 in |25|. 
If we denote by Z' r '(0) the coefficient of x°y r in 

X(Y[a,y]a~ 1 ,x), 
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then 

OC i oo 

tri.n^ ,(o) = tr i^n^ (l)(o) - ( 3 - 13 ) 

i=0 i=0 

In particular if we consider only the coefficient of q$ , i.e. V(ti,T2, ■■■), we get 

tr|^n^ <I,(0) =tr|M (1) naf )(0) , (3.14) 

i=0 i=0 

where M(l) is identified with M(l) ® 1 C V L by means of J£E). 

Remark 3.2 Generalized characters discussed in this section already appeared in the literature. 
in ED ID similar generating functions related to mirror symmetry on the torus resemble (13.711 . 
Also they appear (with finitely many L^(O)'s) in |19|. 

3.3 Bloch-Okounkov's n— point function 

In [3 Bloch and Okounkov introduced another type of generating functions. Let us associate to 
the generalized character V(ti, T2, ...) an n-point function in the following way. First we define 
operators 

. s 1 1 A d y r 

r — l 

Then the (formal) n-point function is defined as 

F(Vl, -,Vn) = V(Ql) ff (yi) ■ ■ ■ v(yn)V(T 1 ,T2, ...)|r 2 =r 3 = ...0. (3.15) 

We mention that m's is a formal variables in contrast to |3] where Ui is a complex variable. 
Proposition 3.1 We have 

r)(qi)cr(yi) ■ ■ ■ cr(y n )V(Ti,T 2 , ■ ■■)\t 2 =t 3 = ...0 = 

Coeff :E otrj^ J7((7i)X(y[V',?/i]^*,xi) •• ■ X(Y[ip,y n ]ip'" ,x n )q^ 0) = 
tr\r a r,( qi HY[^, y^*) ■ ■ ■ o(Y&, y n W)qf (0) , (3.16) 

where L(0) = L(0) - -L. 

Proof: The second equality is trivial since 

X(u, x) = u(wt(w) — n — l)x~ n . 

nez 

For the first equality we prove only in the case n = 1. Let us recall 



V(ti,t 2 ,...) =tr|^ ~[[qf ' 



i=l 

Then 



v(<nW(y)v(T 1 ,T 2 , ...)| T2=T3= ... = 



1 °° r 

trko «(9i)(- + E Z(r) (°)7r)^ 



i(0) 

r! 
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/ 1 00 r 

trta --^C(-r,l/2)^+o( : X(^ e H x)X(^,x) : ) 

\* r=l 



-o(:Jf(^,ar)Jf(^i):) ]gf (0) . (3.17) 



Since the charge operator o(' m X(ip,x)X(ip,x)\) acts as zero on To and 

„»/2 

e* 

r — 1 

the proof follows from 13,51 , 



1 1 y T _ e y/2 



Now let us apply the results from Section 2.4 in the case of vertex operator superalgebra 
T (with the obvious super-modification). Fix m = ip, Vi — ip*, i = l,..,n. Notice that (un- 
like in Section 2.4.) we do not consider an ^-module but rather a subalgebra of T. Since 
X(t/j,tiXi)X(if)* ,Xi) acts on To we can compute the q-trace 

tr|*bX(V>, tixi)X(i>*,xi) ■ ■ ■ X(i>, t n x n )X{ip*, x n )q L{0) , (3.18) 

and a corresponding coefficient of a; , multiplied by T](q): 

G(t u ...,t n ) := (3.19) 
= r)(q)tT\r o(X(i>,tiXi)X(i>*,xi)) ■ ■ ■ o(X(ip, t n x n )X(i>*, x n ))q L(0) . 

Proposition 3.2 Suppose that 13.1<A) converges, in some domain, to a multivalued function 
G{t\, t n ). Then, as complex functions, 

F(Vl, ■■; Vn) = G(tl,..., tn)| t n/2 =e „ w/2 

for every y provided that the right hand side is convergent. 
Proof: Since 

ip(n)tp* = 5 n ,i/ a l, 

for n > j and 

f l/2 



X(iP,tiXi)X(ip\xi) = lX(ip,Uxi)X(i> ,Xi); + 



t-V 



X(Y[iP, Vi }iP*,Xi) = ;Xy,,e vt Xi)X(4>;xi)l + -^—j, 

we can apply Theorem 2.1 and Corollary 2.1 (with minor modifications). The correlation function 
G converges to a multi valued function (cf. Remark 12.21 . However, we can choose a branch 
£»/2 _ e «yi/2 Therefore the Proposition holds. I 

We will need the following trivial (but important) observation. 

Lemma 3.1 Let tp : Vi — > Vi be a vertex operator superalgebra isomorphism and W a subalgebra 
of Vi stable under X(ui,x±) ■■■ X(u„,x„). Then 

tr| w Jr(«i,xi) • ■ ■ X(u n , x n )q Lv i (0) = 

trl^X^K),^!) • • • X(ip(u n ),x n )q Lv 2 {0) . (3.20) 
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Proof: It follows immediately form the property of the trace 

ti'l^w) , ii) ■ • • X(u„, i„)( ( 3" 1 g il '2 (0) = tr\ M X(ui, Xi) ■ ■ ■ X(u„, x n )q Lv ^ {0 \ 

and the fact that 

<pX(ui,a;i)v9 _1 = X((p(ui),xi). 



3.4 q— traces in the case of free fermions 

The boson-fermion correspondence implies (cf. Lemma l.'i.ljl that 

tr|jrX(^, tixi)XW t xi) ■ ■ ■ X(i>, t n x n )X{^*,x n )q L(0) , (3.21) 

is equal to 

tY^X^t^Xia- 1 ^!) ■ ■ ■ X(a,t n x n )X(a-\x n )q L{0 \ (3.22) 
and the corresponding trace for the zero-charge subspace Fa is equal to I3.18H . 

oo 

Now we need some notation. Let (q)oo = — <f) an d 

i>l 

e(t) = ^2q^t n eC[[q 1 / 2 ,t ±1 ]]. 

lit = e 27viv £ C x and q = e 2niT , r € H then 9(t) is the classical theta function. Also we define 
another theta function in the infinite-product form 



hi(r,v) = 



t~ 

i=l i=l 



(9) 2 oc 

This is an entire function of v for fixed r £ EL For a fixed r we will write #ii(n). These theta 
functions can be thought as functions of t (instead of v) but then these are multi valued function 
which we will denote by 8(t) and 0u(f). The following transformation formulas hold: 

e 11 { q t) = -q- 1/2 t- 1 e{t). 

The following result has been known, in some form, for a while. For the completeness we 
include a proof here. 

Theorem 3.1 (a) 

tr\ T X(iP, tiXi)X(i/>*, xi)--- X(i>, t n x n )X(i>*,x n )q Lm (3.23) 
converges to a multi-valued holomorphic function in the domain 

\ti_xi\ > \xi\ > ... > \t„x„\ > \x n \ > \qtixi\ > 0, (3.24) 
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(b) 

r]{q)tr\r X(^, tixi)XW, xi) ■ ■ ■ X(i>, t n x n )X{^*,x n )q L(0) (3.25) 



converges in the domain [3.2$ , an d it has a meromorphic extension to a double cover of 
(C x ) 2 ' 1 equals to 



\\0n (tjXj/tiXi)Ou (xj/xi) 



(3.26) 



Y[6u faxj/xJOn {xj/Uxi) ]J On (U) 

i<Cj i—1 

Proof: Let 

G n (t,x) := X{il),xiti)X{il)*,ti)---X(ip,x n t n )X{<ip*,t n ). 

First notice that, because of the boson-fermion correspondence, we can replace 
X(i(j,xiti)X(i)*,ti) by X(a,x 1 ti)X{a,-\t l ). Also 

tr\ M ( m c[L]G n (t,x)q Lm = (tr\ MW G n (t,x)q LW )(tr c[L] G n (t,x)q Lm ). (3.27) 

We fix a basis for M(l): 

h{-mp ■■■h(~n l ) k >l 



..,ki e N. Al 
h'(-n 1 ) hl ---h'(-ni) h n 



where h = -^=, m > n2 > ... > n; and ki, ki G N. Also we fix the corresponding dual basis 



Now we calculate (as in 28 ) 

tr C [h(-m)]G„(t, x)q L{0) = 

E lAT(h'(-m) k ,G n (t,x)q L ^h{-m) k ) = 

k£N 

Yl{(tjx J /t i x l ) 1/2 - (t :j x j /t i x i y 1/2 )((x :i /x i ) 1/2 - (xj/x t y 1/2 ) 

E 



ITltH 



nc*, i/2 -t-^natjXj/x^ 2 - (w*o~ i/a )(te/*i*o i/2 - {xjiuxi)- 1 '*) 

i—1 



^(ft'(-ro)*:G»(t,x):ft(-m)*>, 



klm k 



Then 



E^^'(-™) fc : G »(*>*): ft (-™) fe > = 

fc€N 

fcGN ' \ r=l 

exp > (avir) + £ r /if — m) ) = 
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E fci^r E |(-^ E^') m + ^ m )>'(-™) fc ,ex P ( - 

fc£N i=0 r \ 

_ \ 

^(a; s t s ) m + x s m J a(myq mk h(— m) k ) — 
•=i / 

k\(h'(-m) k , h(-m) k ) (— £ E r , s Ifrffr + 



oo oo 



EE ,„/./,.! (/,•-/)!/!* 

i=0 fc=i v ' 



mk 



E 

7=0 



/ ]_ y t ^ x > Xr_ _ t™x™ _ x r \i 



m £~tr,s tV^X 1 ? 1 x s x s (^i™ / (1 — q m ) 1 

i!(l - g m ) 



exp I + _!_. (3.28) 

l t^x™ x a x a t^x™ l-q m 1 - q m 



r.s 



So far we calculated the trace on the space C[/i(— m)]. Since 

M(l) S5 ®~ =1 C[/i(-m)], 



then 



Therefore 



tr| M( i)A = JT trj c[M _ m)] A 



n 



1 4-m m fm m n 
± ,v r a/ r L s a> s , -i \ y 

m=l \ ■ s s r 



TT exp --(i^ + ^^ + l)-^ = (3.29) 

11 \ m^t^x™ t™x™ 'l-q m l 



for 1 < r < n and 1 < s < n. Thus 

n 8<J ((w^) 1/2 - (w^r 1/2 xfe/^) 1/2 - fe/xQ- 1 / 2 ) 
nr=i& 1/2 - ^ 1/2 )n i<3 (fe^M) i/2 - (w*o- i/2 )((*j/*«zo i/a - fe/te)- i/2 ) 



n e H~™E(; 



. m ^ tl^x™ x s x s ti^x™ 1 - g m / 1 - g" 

m— 1 \ r.s ' 



Ili<3 Ou(tjXj/tiXi)Oll(xj/Xi) ^ ^ 



n«i fu (tjXj/xiWu {xj/uxi) nr =1 #n (*<)(«) 

If we assume that 



> |a;i| > ... > \t n x n \ > \x n \ > \qhxi] > 0, 

then the infinite product 13.291 converges uniformly to a holomorphic function inside (13. 241 . This 
is a proof of (b). For (a) we have to calculate (cf. 13.271 1 ) 



tr\ c[ L]G n {t,x)q Lm ='£t?-.-Cq S t 



n£2 



which is convergent for \q\ < 1. 
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From Theorem 13. II it follows that 

n(q)tr\ To X{^,t 1 x 1 )X{^ ,xi) ■ ■ ■ X(i/>,t n x n )X(il>* ,x n )q LW 

has a meromorphic continuation (after we specify branches) to (C x ) 2n with the set of singularities 
determined by zeros of 

n 

Y{e 11 {U)Y{e 11 (t J x ] /x 1 )9ii(x 3 /Ux l ). (3.31) 

i — l i<3 

Let x — (xi, x n ). Denote by A n (x) the set of all n-tuples (t\, ■■■,t n ) such that (ti, t n , xi, x n ) 
satisfies inequality H3.24|l . 

3.5 Another description of the q— difference equations 

From Theorem 13. lf b) it follows: 

Corollary 3.1 Let \q\ < 1. Then 

tr\ ro r](q)X(a, tx)X(a~ 1 , x)q L(0) (3.32) 

converges in the domain tjt > \t\ > 1. Moreover, has a meromorphic extension to a double 

cover of C x equal to 

1 



0il(i)' 

Corollary 13. II is essentially Theorem 6.5 in |[|]. 
Remark 3.3 According to 0, the following formula holds 



Hvi,-,Vn)= £ — 7TT f — VTT f 1 'Y' 1 ' (3 - 33) 

where ti = e Vi . Let us denote by G(t%, t„) the coefficient of x° in the expansion of I3.25H . 
Because of Proposition 13.21 

Tijji, ...,y n ) = G{ti, ...,t n )\ t n/2 =eny . /2 . 

If we switch to complex variables a Laurent expansion of 13.26H (inside a certain domain) gives 
us an integral representation of G(ti, t n ). Let Ai be an annulus inside the z—plane such that 
Ai is contained in (13.2411 . and d is a smooth curve inside Ai. Then the Cauchy's integral formula 
for several complex variables implies 



Y\ 01 1 (tj X j /ti Xi ) 0! l (x j /. 

i<3 

n 

Y\6 1 i{t ] x J /x 1 )6 1 i{x ] /t i x 1 )Y\6 11 {U) 

i<.j i—l 

This integral is hard to calculate explicitly. 
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Let us consider the case n = 2 in more details. Let x = — . Assume that g 



n /lit tit 



r e 



and | £2 j > tij. Then the set of singularities 13.311 is ordered (with respect to the absolute value) 
in the following way: 



> 







q k+1 




q k + 2 




> |*i I > 




> \qti\ > 




t 2 


t 2 


t 2 



for some k £ Z. Suppose for simplicity that k = 0, i.e 

1 

*2 



> 



> *1 > 



> 1^1 I > 



> 



(3.35) 



(3.36) 



S- 

»2 



Put n = \qti\, r 2 
radius r. From Remark 

Corollary 3.2 Let 



and r3 = \t%\. Let us denote by C r [z) G C a circle with center at z and 
-{.31 in the case when n — 2, it follows 



N > 



t 2 



Then 



where 



G(t u t 2 ) = 



dx 9 11 (t2x/t 1 )6n(x) 
(ti)9u(t 2 ) Jc (q) 27ria; 0ii(<2a;)0ii(a;/ii) ! 



ai 



G(qti,t 2 ) = 



where 



Since 



«_> 



G(qti,t 2 ) = 



G(qti,t2,x) — —q 1 ' 2 tit2G(ti,t2,x), 

dx Qii{t2x/t\)9\i(x) 



q ll2 Ut 2 



9u(t 1 )9 11 (t 2 ) /c Q (0) 27ria;6lii(t2a;)6lii(a:/ti) : 



(3.37) 
(3.38) 



Let us try to determine the behavior of 13.381 under the elliptic transformation t\ i— > qt\ . Now, 

1 [ dx 9n(t 2 x/qti)9n(x) 

hi(qti)6n(t 2 ) J Ca {0) 2-Kix 9n(t 2 x)9 11 (x/qt 1 ) ' 



Now from 13.361 and Cauchy's theorem we have 

l 1/2 t!t 2 I dx 8ii(t 2 x/ti)9 n (x) 



G(qh,t 2 ) 



-q 1/2 tit 2 



hi(ti)9u(t 2 ) Jc a2{0) 2 ™x 9u(t 2 x)9n{x/t 1 ) 
1 f dx 8n(t 2 x/ti)9ii(x 



9u(ti)9u(t 2 ) J Cai (o) 2 ™x 9 11 (t 2 x)9 11 {x/t 1 ) 
dx 9u(t 2 x/ti)8ii(x) 



+ 



9ii(ti)9u(t 2 ) Jce(J-) 2mx9 11 {t 2 x)9 11 (x/t 1 ) 

(3.39) 

for some small e > (cf. Figure where dotted circles have radius ai and a 2 and the small 
circle is centered at Let us evaluate the second integral on the right hand side of 1)3.39(1 . 
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Figure 1: Contours of integration 



Lemma 3.2 For small e > 
(a) 



dx 8ii(t2x/ti)8n(x) 



Mti)0ii(i 2 ) Tc t ( qtl ) 2™ 9 11 {t 2 x)9 11 (x/t 1 ) 0n(tit 2 ) ' 

and 

(b) 

dx 9 11 (t 2 x/t 1 )9 11 {x) _ 1 



0ii(ii)6>n(fa) J c ,f jli 27ria; M^M^Ai) 9 u (t 1 t 2 ) ' 



(3.40) 



(3.41) 



Proof: The proof of (b) is trivial once we prove (a). Note first that the two functions in l|3.4()^ 
are multivalued. Since we can extract on both sides a term s/t\ti we may think of I3.4UH as an 
equality of single-valued functions. 

dx 0i\{t2x/t\)6\i{x) _ 



9n(ii)0n(fe) Jc 3 2TTix9 11 (t 2 x)8 11 (x/t 1 ) 
1 



0n(ii)0ii(fe) 

dx 9ii{t2x/ti)9ii(x)(q)'j _ 

c 3 2mx(l - 9ll (t a x) IW1 - ^)(1 ^ ^)((*A0 1/2 - (x/t,)-^) ~ 

gll(g*2)flll(gtl)(g)gc = 

9ii(ti)9ii(t2)9n(qt 1 t2)U i >i(l ~ - 1 l )(l 1/2 - 1~ 1/2 ) 

' (3.42) 



?ii(ta)0ii(ti)(ff)2o 1 



6>ii (<i)(9ii(t2)eii(tit2)(3)io 0n(iife) ' 

I 

Lemma |3.2I and 13.3911 imply that the g-difference equation for G(ti,t 2 ) depends on the 
one-point function G{t\t 2 ), i.e. 

G(qt!,t 2 ) = -q 1/2 t 1 t2(G(t 1 ,t 2 ) - G(tit 2 )). 
In the next section we will prove the same formula by using purely algebraic methods. 
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3.6 q— difference equations for Bloch-Okounkov n— point functions 

In this part we analyze the difference equations for G(ii, t n ) (at least when n = 2) which gives 
us elliptic transformation formulas for the multi valued function G(ti, ...,t n ), i.e. transformation 
formulas for G(qti, t n ). If we put 



J\0u (tjXj/tiXi)9u (xj/a 



H it i j ...j tn •) x i j . . . j ) 



then we have 

H(qti, t n ,xi, x n ) = -q 1/2 t 1 t 2 ...t n H(t 1 , t„, xi, x„). (3.43) 
where we used the formula 

e 11 (qt) = - q - 1/2 t- 1 e 11 (t). 

We are aiming for a similar formula as in the case of G(ti, ...,t n ). For simplicity let us consider 
first the 1 and 2-point functions. 

We start with the 1-point function 

G(h) = r?(g)tr|^ X{$, Uxi)X(i>* , xi)q t{0) . (3.44) 

Then we have 
Proposition 3.3 

G(gti) + q 1/2 t 1 G{t 1 ) = 6 1/2 (qh). (3.45) 

Proof: First note that 

[X(ip,t 1 xi),X(iP*,xi)] = S l/2 (ti). 

Also, we introduce a "control" variable y a C[y, y _1 ]-valued operator y H , where H(0) is the 
charge operator introduced earlier. 

r,(q)tr\fX(4,, t lXl )X(r,xi)y HW q m) = 

»?(ff)trMi/ a (ti)y* ( V <0) - V (q)tr\ T X(r,x l )X(i>,t,x 1 )y H ^q m = 
r/(?)trM 1/2 (*i)5 £(0) -tr|^WVi)y K( ° H V (0) ^ — ) = 

^trkJi/afa)^ - ^(q)tr|^(^, — )X(V*, xi)y H(0) - V (0) - (3.46) 

Now if we extract the coefficient of t/° in (13.4611 we obtain 

G(ti) = *i/a(ti)-3 _1/3 tiG(ti/?). 
After multiplying with q D *i we obtain the desired formula. I 

Note that l|3.45[l is a formal analogue of the difference equation for one-point functions ob- 
tained in P). By using the simple identity (1 — t)5i/ 2 {t) = 0. we obtain 

(1 - qtx){G{qtx) + g 1/2 tiG(ti)) = 0. 
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Note that in the equation above we cannot remove the term (1 — qt\) because this is a formal 
expression. On the other hand this equation gives us meromorphic continuation on (a double 
cover of) C x (provided that G(ti) is meromorphic in Q, and such that q z ■ Q — C x ). Now we 
consider transition from the formal variables to the complex numbers. From 113.451 it follows 

G(qti) + q 1/2 tiG(h) = (q Dt i + g 1/2 ^)G(£i) = 
= (l + g Ct ig- 1/2 irV /2 iiG(ti) = 

= S 1/2 (qt 1 ). (3.47) 
Now we formally invert the operator (1 + q Dtl 9 _ ijf 1 ) by using a formal expansion 



— ^ = Y(~l) n A n . 
1+A ^ y ' 



Then we have 



ql/2tlG{tl) = Tw^HT Mqtl) = 

' S 1/2 {qti). (3-48) 



Therefore 



1 + q-VH-iq 13 *! 



G(ti) = q~ 1/2 t^ (-ir(?- 3/2 tr 1 «^) m 5i/2(«ti; 

m>0 



-1/2,-1 I -3/2,-1 D t . \m n.n 

1 1 h 2^ 2^,{l ' h 1 H ) 1 *i 

neZ+l/2 m>0 
l(m+l) 



E E (-l) m Q nl ^q- kim+1) h k - 1/2 , (3.49) 



fe£Zm>0 

where n~ m — — fc + 1/2. Now, since 



m(m + l) (m-fc + 1/2) 2 (fc + 1/2) 2 

~ /C( 771 "-f* J- ) — — — — . 



we obtain 



g(«i) = E <r 7 « — ^~ E c^r* — 5 — • ( 3 - 5 °) 

feez m>0 
By using the symmetry on the right hand side of 13.5011 and after summing the geometric series 
we obtain: 

ceo - £ i-vrf** + ■ + J^r) . <«.) 

m>0 v 1 7 

If we switch to complex variables then the Note that H3.51H is a convergent series when 1 < \t\ < 
jij , and it can be analytically extended (by using j3.5ip ) to all values t 7^ q n , n € Z. Thus 

G(*0 = *r 1/2 E (-l) m <T^ J— i . (3.52) 

We claim that G(ti) in 13.521 1 is equal to 13.441 (remember we just find a solution of 13.3H '). For 
this apply Lemma 18 . 21 from Appendix A. 
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Now, if we combine formula 13.521 with Corollary 13. II we see that the reciprocal theta 
function can be written as an infinite sum (with the same set of poles). The similar formula 
follows from the denominator formula for N — 2 superconformal algebra. 

Now we discuss two point functions. Again the aim is to derive q-difference equation for the 
two point function by using a formal variable approach. 

Proposition 3.4 We have 

G{qtiM) + q 1/2 tit2G{t 1 ,t 2 ) = S 1/2 {qt 1 )G{t 2 ) + g 1/2 ti£ 2 G(M 2 ). (3.53) 

Moreover we have a locality equation: 

(1 - gii* 2 )(l - qti) {G(qti,t 2 ) + q 1/2 ht 2 G(ti,t 2 ) - q 1/2 t 1 t 2 G{t 1 t 2 )] = 0. (3.54) 

Proof: The proof is very similar as the one in Proposition 13.31 . So we omit details. By using 
the Jacobi identity 

tT\ r ri(q)o(X(i},tix 1 )X(^* ,xi))o(X(ip,t 2 x 2 )X(ip* ,x 2 ))y H{0) q L{0) + 
tr\rv(q)o(Xy>, — W, x 1 ))o(X^,t 2 x 2 )X(r, x 2 ))y HW ^ q L{0) = 

q 

tT\rviq)Si/2iti)X^,t 2 x 2 )X(i>*,x 2 )y Hm q Lm - 

tr\ T X{^Mt2X 1 )X{^\x 1 )y H ^\ Uo) +tr\ F 8 V2 {t 1 t 2 )y H{0) q L(0 K (3.55) 
If we extract coefficient of (xix 2 )° in 13.551 . by using Theorem 3.1, we obtain 

-V2, . n ( *i 



Si/ 2 (ti)G(t 2 ) - G(ht 2 ) + S 1/2 {tit a ). (3.56) 

If we act by q Dt we obtain 

G(qt!,t 2 ) + q 1/2 t 1 t 2 G(t 1 ,t 2 ) = 

(7 1/2 tit2G(tit2) - 5 1/2 (qt 1 t 2 ) + 5 1/2 (qti)G(t 2 ) + S 1/2 (qt 1 t 2 ) = 

g 1/2 fif 2 G(tit2) + 5 1/2 {qti)G(t 2 ). (3.57) 



By using this method we can derive similar formulas for an arbitrary n-point function. A 
similar theorem has been proven in [3] by using different methods. 



Theorem 3.2 Formally 
(a) 



G{t u ...,t n ) + q~ 1/2 ^11**) G^,i 2 ,...,t„ 

n-1 

~ /, z2 G(t 1 t il ■ ■ -ti 3 , . . . ,i ix , . . . ,i ik , . . .) + 

8=1 Kti 1 <...<ti s <n 

+S 1/2 (U 1 ■ ■ ■ U,)G{- iu ■■■)■ + 8 1/2 (ti)G(t 2 , t n ). (3.58) 
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(b) 



Yl (l-qtiti 1 ...t is )(G(qti,...,t n )+ (3.59) 

l<il < . . .z s <n \ 

71-1 \ 

i 1/2 Ii u T, E (-i) s G(t 1 t il ---t is ,...i il ,...,t is ,...) )=o. 



0<s<n — l,l<ii<...z s <n 
n n—X 



i=l s-0 Kii<...<i 3 <Ti 

(cj Lei |g| < 1. 27ien G(t\, ...,t n ) converges inside an open subset 



n 

O«,on{(ti,..,i») : | J] ti| < -}, 



and it /ias a meromorphic continuation on a double cover of (C x ) n tuit/i t/ie set of poles 
contained in {(ti, t n ) : <? m tn ■■ - ti„ = 1, 1 < ii < ... < i s < n, m £ Z}. 

Proof: The idea is essentially the same as in (13.5511 . We prove (a) directly and use the induction 
for (b). Actually (b) follows if 

n 

G(qh, ...,t n ) + q 1/2 J] UG(ti,. ..,*„) = 

n 71—1 

5 1/2 n^E E (-i)- 1 G(t l *i 1 ---ti.,...,fc 1> ---,^...o 

i— 1 s — 1 l<i 1 <, ,.<i s <n 

+ E E a « jfc^ii-'-^JGO..,^,...^,...), (3.60) 

where ctj 1 ,...,j k are some constant. By multiplying 13.601 by 

Yl (1-qtiU, ■■■ti B ) 



0<s<n-l 
Kil<-..<i s <n 



we get 

For n — 1 and n = 2 I3.60H holds. Suppose that it holds for every k < n — 1. Then 



fc 



E E a s G{qt\t ix ■■•ti s ,...,ii 1 ,...,ti s ,...) = 



3=1 Ki 1 <...<i„<fc 
A' 



Et- 1 )'^ 1 E E L )(-ir« m G(iit il --.i ir ,... I t il) ...t ir) ...) + 

r— 1 Kji<...<j r <fe m— 1 \ / 

(3.61) 

Where dots stand for terms that involve 5{qti 1 • • • ti k ). In particular if for every i, at = 1 then 



E(;)(-ir-=-i- 
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and 



G(qtiU 1 ■ ■ ■ U a , . . . . . . ,U a , . . .) 



s — \ ii <...<i s <fe 
k 

-£(- 1 ) r ~ 1 E G(t 1 t il ■ ■ - u r 

r=l Kj 1 <...<j r <k 

Let o stands for the operator Coff, z .o.... I ,o • 



• , ti 1 , • ■ ■ , £j r ...)+.. 



(3.62) 



Therefore 



tr 



tr|jr o(X(^,x 1 t 1 )X(V'*,ti)---X(^, a;n t„)X(V*,t„))g £(0) y H(0) = 
tr| r o([X{i>, Xiti),X(1>* , xi) ■ ■ ■ X(i>, x n t n )X{^* , x n )])q L{0) y H{0) 
tr| r o(X(V*,xi)---X(t/ ) ,a; n t n )X(V*,a : „)X(V,siti))g L(0) y J/(0) . 

r o(X{ip, Xit^Xty^Xi) ■ ■ ■ X(i>, t n x n )X(^,x n ))q m y H(0) + 
T o(X(rP, El*L)x{iP*,xi) ■ ■ ■ X(i>,x n t n )X(i>* ,x n ))q m y H(0) - X = 
T o(S 1/2 (ti)X(iP, t 2 x 2 )X(r ,x 2 )... X{^, t n x n )X(i,*, x n ))q Lm y H 



(3.63) 



(0) 



+ 



:\f o(X(tp* ,x 1 )X(^,t 2 x 2 )S 1/2 

■\MX(r,tl) • ■ • X(l/>, t n x n )S 1/2 (^-) )q Lm y H(0) . 
Now in 13.6411 we use a formal delta function substitution property (cf. [Ill) 



t\x\ 

x 2 

t\X 

X 



■ ■ ■ X(^,t n x n )X(r,Xn))q L(0) y Hm + 



X(u, t 2 x 2 )5 1/2 ( -!— - \ = X(u, t^x^S ( — - 



x 2 



t 1 X 1 



■r 2 



(3.64) 



(3.65) 



This fact was used in Part I |25| as well. Then 

tr|^ o{X(i;,t 1 x 1 )X(i>'" ,xi) ■ ■ ■ X{^,t n x n )X(^* ,x n ))q LW y Hm + 
tt\r o(X(^,— W • ■ X(i,,t n x n )X(r ,x„))g £(0) y H(0) - 1 = 

tr|y o(*i /a (ti)X(V, t 2 a;2)X(^*, x 2 ) ■ ■ ■ X(V», t n x n )X(i/> m , x n ))q L{0) y H{0) + 

tr|y o(X(^',a!i)Jf(V,titaa:i) • ■ • t„^„)X(^*, x„))g i(0) y H(0) + . . . 

tr|jr o(X(^,o;i)X(V,t 2 a; 2 )---X(V,ti*„a;i))g £(0) y H(0) . (3.66) 

We move the term X[%j), titiXi), i = 2, ...,n to the right hand side of 13.6611 to the left by 
anticommuting them with operators X(ip* ,Xj), j < i. This anticommuting produces more terms, 
etc. If we repeat procedure of anticommuting and delta function substitution for each term 
(except the first term) on the right hand side of (13.661 we obtain 

tr\ T o(X(i>* ,xi)xy>,titaxi) ■ ■ ■ X(^,t n x n )X(^ ,x n ))q L{0) y H{0) = 
\x\t S 1/2 (t 1 t 2 )o(- ■ ■ X(i>,t n X n )X(i>* ,x n ))q Lm y HW - 

tr ^ o(X(^,t 1 t 2 x 1 )X(^* , Xl ) ■ ■ ■ X{i>,t n x n )X{i>* ,x n ))q Lm y Hm , (3.67) 
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for the second term, 

tr|jr o{X{ll>\ Xi)X{^, t 2 x 2 ) ■ ■ ■ X(i>, tiUzi) ■ ■ ■ X{i>, t n x n )X{^*,x n ))q L{0) y H(0) = 

i-1 

-J2 J2 tr l^ o(X(1),t ai •••t Sr tiXi)X(ip*,xi)--- 

r=l l<t,j <...<t Sr <i 

■■.X(iP,t nXn )X(iP*,x n ))q L(0) y Hm 
+tr|jr 8 1 / 2 (t 1 ■ ■ ■ i si ■ ■ ■ is T ■ ■ ■ U) ■ 

o(X(i>,t sl x sl )X(i>*,x sl ) ■ ■ ■ X(^,t Sr x Sr ) ■ ■ ■ X(^,x Sr ))q Lm y H(0) (3.68) 
for the i + 1-th term, etc. Combining 13.661 , 13.671 and 13.681 we obtain 

n t 

G(t u ...,t n )+q- 1,2 T\ UG(-,t 2 ,...,t n ) = 

i>l H 

n-1 

— 2^ G{t\ti 1 • • • ti a , . . . , ti t , . . . , ii h , . . .) + 

s = l Ktjj <...<t; s 

+5 l/2 (ti 1 ---ti.)G(...,i il ,...ii....). (3.69) 
This is exactly 13.581 . Now if we act by q *i on 13.691 we obtain 

G(qti, ...,t n ) + q 1/2 h . ..t n G(t!,t 2 , ...,*„) = 

-^2 Giqtit^ ...t is , ,ti k ,"-) + 

s>l Ktjj <...<t is 

+%t il ...t ij> )G(...£s 1 ... &,...)• (3.70) 
If we apply now 13.621 and the inductive assumption we obtain 

G(qti,...,t n ) + q 1/2 ti...t n G(ti,t 2 ,...,t n ) = 
l 

'1 

i — 1 s>l l<ii<. ..<i s <n 

+5(qt h ...ti,)G(...i H ...ii a ...) + ... (3.71) 



1,1/2 II ^X] H (-i) s ~ 1 G(txu 1 ...t is ,...i il ...t i 



which is 13,601 , There are essentially two ways of proving (c). First which uses an explicit 
g-expansion for G(ti, ...,t n ) and the second that uses Theorem 13.11 The first method is much 
harder though. Let us explain the second proof. The crucial observation is that 

T„nfi„, cA„(i) (3.72) 

holds for every x where 

n 

T„ = {(ti,...,t„) : | JJti| < -}. 

i=i ^ 

This is easy to prove by induction. Since the 2n-point function converges absolutely inside fi n ,™ 
and it has an expansion in terms of x^s, all terms converge absolutely, hence the zeroth term 
G(ti, ...,tn) converges absolutely inside T n PI Zl ni o. Analytic extension and information about 
poles follows from the part (b) by using the induction. 
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Now let us try to solve 113.5611 for n — 2. For n > 3 it is possible to obtain similar explicit 
formula. Because, the expansion in ti's converges in a tiny domain, i.e., Y„nf2 n ,o, these formulas 
are hard to express in a nice way (meaning, that the poles and domain of convergence are visible) . 
On the contrary, product formulas have very nice form (cf. 3;). 

Proposition 3.5 (a) 13.561 ) has the unique solution in C[[q, if 1 , t J 1 ]] that satisfy G(fa,fa) = 
G{t 2 ,t\) and it is given by 



G(t 1 ,t 2 ) 



V2 VL!^ 1 ^ 1 ~1 



m>0 

-l.-k 



(1 - (iifa)-^™)^ - (fafa)9 m+1 ) 



1/2 \ "> tl ^2 

q \l-k\ 



(^r 1/2 E i 

(k,l)£Z 2 ,k^l 



I X, (-1)^ 5 2 — 9 (~ 1) 9 

\m=0 m=0 / 

(b) Let \q\ < 1. 27ien G{ti,t2) absolutely converges inside 

^2,0 n {(ti,t 2 ) : |*x*2 1 < A}. 

Moreover, it has a meromorphic extension to a double cover o/(C x ) 2 with the set of poles 
contained in 

(C x ) 2 \{(ti,t 2 ) : nt 2 q m = l,t iq m = l,t 2 q m = l,m G Z}. 

Proof: The fact about the uniqueness follows from Lemma 18.41 from the Appendix A. Proof of 
part (a) is quite lengthy, but straightforward and essentially uses the same procedure as in the 
case of 1-point function. Here is the main idea. Formula 113.5611 can be written as 



g 1/2 tit 2 G(ti, fa) = ri 5- (<5 1/2 (gti)G(fa) + g 1/2 tifaG(tifa)) . (3.73) 

The fraction in the above formula can be expanded as in the n = 1 case. The rest is manipulations 
with g-series. 

We can prove part (b) in two different ways. The fact about the convergence of G(ti,fa) 
follows either from the previous theorem (which holds for every n) or by using explicit expression 
i3~73t. Hence by using this fact we can meromorphically extend 13.561 to remaining values of t 
by using g-difference for G(qti, fa) and G(q~ 1 ti,t 2 ) and then by using G(ti,qt 2 ) and G{t\, q~ x t 2 ) 
(which is obtained from 13.561 because G(ti,fa) = G(fa,ti)). From the g-difference equation it 
follows that we can not extend the function to those values for which either G(tifa) has poles or 
G(t\) and G(t 2 ) hence the statement follows. I 

It is worth mentioning that (as in the n = 1 case) we encounter certain g-series that have 
singular parts consisting of single terms in the g-expansion. Once we combined several of these 
series all singular terms cancel. This should be compared with some results from Part I |25|. Let 
us illustrate this using a particular example. An incomplete ^-function is an expression of the 
form 

(m + l-2fc)(m + 2) 

Mi) = (- 1 ) m 9 

m=0 
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for some k 6 Z. It is clear that Afc(g) £ Z[[g]] for fe < 0. Also it is not hard to show that 

1 i , i 1 m(m + 2fc+l) 1 

Mi) = i + E + 9 — 5 — G i + z ^ 

" m=0 * 

for fc > 0. Then 

(g) - g 1 ' fc| A max(feji) (g) G Z[[g]], 
for every choice of constants l,k G Z, k I. This fact is used in proving 13,731 . Also notice that 
coefficients of f^ fc ~ 1,/2 t^ !_1 ' /2 in 13.731 are given by 

A klt = J2 (-ir +n 9 ln+m+ T +m+2> e Z [[ 9 ]], 

m,n>0 

for every k,l € Z, k ^ I. Some related g-series were studied in |15|. 

3.7 General case 

So far we only discussed charge subspace, i.e. J~o C T . However T m is also a ©-module. 
In vertex operator algebra language, T m is an irreducible .7-o-module (M(l) = J-"o — £i cf. 
0). Let us denote by H m (ti, t n ,x 1 , i„) the restriction of 13.231 to the space T m and with 
Gm{ti, t n ) the corresponding Bloch-Okounkov n-point function. For every m G Z we have 

H m (qti, ...,t n ,xi,..., x„) = — g' I+1/2 iit2 • ■ • t n H m (ti, t n ,xi,...,x n ). 

Also, 

C? m (tl, tn) = q" 2 ^ ■ ■ ■ &G{t U ...,t»). 
Hence the "total" Bloch-Okounkov n-point function obtained by taking the trace over IF, and 
it is given by 

/ n \ m m 

E? m2/2 (II*i) G(t u ...,t n ) = 9(j[[u)G(jt 1 ,...,t n ). 

4 g— traces for D ; Majorana fermion case 
4.1 Iterated 2n— point functions 

Let us recall some notation from first part 125) . Denote by F the Fock fermionic space associated 
to a free (or Majorana) fermion. Let <p = y>(— 1/2)1 G F, such that 



X(<p,x) = (^„a; ™ x/ 



1/2 

and 



[95m, 'fin] = <5m + n,0- 

The aim is to study g-graded traces of the form: 

tr\ r X(Y[<p, yi ]<p, xi) • • • x„)g £(0) . (4.1) 

If we modify Corollary 12 . 21 and Remark |2,2l such that it applies in the vertex operator superalge- 
bra case it follows that all the information about 114.11 can be obtained by studying the following 
2n-point function: 

ti\ F X((p, tiXi)X((p, xi) ■ ■ ■ X(ifi, t n x n )X(<p, x n )q L{0) . (4.2) 



30 



Again, once we switch to the complex variables 14.21 is multi-valued with respect to ti variables. 

It is well-known that one can obtain a free fermion from a pair of charged fermions. Let 
£( n ) = <Kn)+|* W ) for n e z + i _ Then 

[v?(m), (p(n)] = <5 m +„,o, 

hence computing 14.11 . up to the character, reduces to the case of correlation functions with 
charged fermions. Hence, it is not hard to see that the correlation function l|4.2[l converges inside 

|*ixi| > Ixil > ... > |x n | > \qtiX\\ > 0, 



to a (multivalued) analytic function. Moreover, l|4.2|l has a meromorphic continuation to a double 
cover of (C x ) 2n . Another way of proving convergence, is by explicit calculations. 
We present the calculations for n — 2. 

tr\ F X((p, tiXi)X(<p, Xi)X(ip, t 2 x 2 )X(ip, x 2 )q L ^ + 
tr\ F X((p, tiXi)X(ip, xi) ■ ■ ■ X(tp, t n x n )X((p, x n )q L ^ = 
ti\ F [X(ip,tiXi),X(tp, xi) ■ • ■ X(if,t n x n )X(ip, x n )]q L(0) = 

tr^i/2 f~) X(^,t 2 x 2 )X(^,x 2 )q L(0) -tr\ F S 1/2 (j^j X(<p,Xi)X(<p,x 2 )q m 



+tv\ F S 1/2 X{vM)X{ V ,t 2 x 2 )q um ■ (4.3) 

tlXl ' 



Hence we get 



tr| F A(<^, t 1 X 1 )X(ip, xi) ■ ■ ■ X(tp, t n X n )X(ip, x n )q L{0) = 
G(ti)G(t 2 )tr| F g £(0) - G{t 1 x 1 /t 2 x 2 )G{x 1 /x 2 )tr\ F q L(0) 
+G(t 1 x 1 /x 2 )G(t 1 /t 2 x 2 )tT F q L(0 \ (4.4) 

where G(t) is defined as in l|3~50jl and 



tr\ F q 



V(qf 



v(q 2 )ri(<i 1/2 ) ' 

As before we consider and the corresponding Bloch-Okounkov n-point function 

D(ti,...,t„) := (4.5) 

= ^^ 2 y <?1/2 ) t r| F o(A(y, tiXl)X( V , xj)) • ■ • o(X(v», tnX n )X{tp, x n ))q L(0) . 

4.2 q— difference equations for n— point functions 

Here we obtain q-difference equations for 14.51 . 
Theorem 4.1 (a) 

D(t!,...,t„) = -D(—,t 2 ,...,t n ) + S 1/2 (ti)D(t 2 , ...,t n ) (4.6) 

n-l 

EE E D<w£...t2 

»=1 Kil<...<ia<n €i ,...,€ is £{ — 1,1} 

+Si/ 2 (t il ■ ■ ■ *i 3 )-D(— i is ,...). 
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(b) 



, . . . ,i a :e,- 



J] (l-qt% ...t%)[D(qt 1 ,...,t n )+ (4.7) 

,..., E ,,e{-i,i} V 



EE E sfttf-C 

s=o i<i 1 ,...,i,<iK, 1 f ia e{-i,i) 

(cj D{t\, i n ) converges inside an open subset ofQ„,o, and it has a meromorphic continuation 
on a double cover o/(C x ) n with the set of poles 

{tH ■••t? m = 1,1 <ii < ••• <in < n S Z>o,m <= Z, ej - G {-1,1}}. 

Proof: Part (a) is essentially the same as the proof of Theorem 13.21 The only difference are 
factors a G {1,-1}. Part (b) follows directly from (a) by acting with q Dtl . Part (c) follows from 
(a) and a similar argument as in Theorem 13.21 



5 q— traces for V + ; bosonic case 
5.1 1— and 2— point functions 

For k > 2, let us define a normalized fe-th Eisenstein series 

— B 2 °° 

G 2k (q) = -j^y- + (27— I)l E ^-iWg" 6 C[[q]], 

where a k is the sum of the fe-powers of the divisors of n and i?2fc a Bernoulli number. If q — e 2nl,T 
we will write G2fc(i~) instead of G2fc(g). Also we define the meromorphic functions 



and 



for k > 2. Hence 



. 1 v-^ 1 1 z 

C(z,t) = -+ ^ + - + — . 

z * — ' z — UJ CO LO z 

oJgZ+Zr 

. , 1 1 1 

P2(Z, T) = - 5 + > 7 j 

u6Z+Zt v ' 



Pk+i{z,r) = ~^7^P fc ( 2 '' r )' 



d 



fc-2 



(-l) k (k-l)\p k (z,r) = f^-J p 2 (z,r). 

((z,t) is the Weierstrass zeta-function, with the property: 

r 



C(2 + r,r) = C(2,r) + 2C^-,r 

For k > 2, pk is elliptic functions ; p2(z, t) is usually called the Weierstrass p-function. In the 
annulus 

< \z\ < min m , n |m + nr|, 
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pk has a Laurent expansion 

1 . 4-l\ 

| G 2 n+2("r)z 



p*(^)=^+(-D-f:(?_ + 1 : 

n=l V 



2n + 1 \ ^ : \ 2n+2-fc 



where G2n+2(i~) are genuine Eisenstein series (not normalized). In parallel with G^kil) it is 
convenient to introduce 



n=l V 



G 2 „+2(g)z + '~ fe . (5.1) 



Then, for fc > 2, we have 

(2TTi) k pk(2-Kiz,T) = p k (z,r). (5.2) 

The Weierstrass ("-function and pk functions are closely related to certain Jacobi modular forms. 
Let us fix (formal) Pfe+i-series for k > 1 

^(..^(g^+ '-'CT'' )- <w> 

Then (15.311 converges inside 1 > \t\ > \q\ if i and q are complex variables. If we let t — e 2n ' v , 
g = e 2 ™ 7 " then in the limit (cf. 

P 2(e 2 ^,a) = 7 ? ^(p2( y ,r) + G2(r)), 

(27Tl)^ 

(Z7TI)' 

for i > 3. Note that we cannot derive relation 15.41 purely by using formal variables. 
The aim is to find the transformation properties of 

- — tr\ M X(u 1 ,t 1 xi)X(vi,x 1 )X(u 2 ,t 2 x 2 )X(v 2 , x 2 )q L{0) (5.5) 
ch M (g) 

under the elliptic transformation: 

ti i— > qti and t 2 t— * qt 2 . 

(because of the symmetry we consider only t 1 t—> gti). 

We will assume that (12.311 holds, i.e., that for every i 7^ j and n > 0: 

Ui(n)vj = c tlii „. j <5 wt ( tl .)_|_ wt („_.)_i in l 

Note that this condition implies 

for every n > 0, where c UyV , n £ C. 
Remark 5.1 

u(n)v = ci5 wt ( u ) +wt (^)_ 1:W l (5.6) 

and 

u[n]v = C25 wt ( u ) + wt(u)-i,nl (5.7) 
rarely hold simultaneously because of the different grading. More precisely, if wt(lt) and wt(w) > 1 
then both 15.611 and 15. 7> hold if and only if wt(u) = wt(w) = 1. This can be shown by writing 
explicitly v[n] in terms of v(i)'s (and vice-versa). 
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Now suppose that wt(w),wt(w) > 1 and that for every i 

o(ui)\ M = 0, (5.8) 

The iterated 2-point function and its zeroth term is easy to compute. From the Jacobi identity 
for X -operators (cf. [23] and |25|) 

|xh„,,^. M .w(^i) Wh *.«) 

it follows 

tr\MX(u, tixi)X(v, = 

tr\ M X(u, t -^)X(v,x 1 )q Lm + Y, S ^ 1 ^\MD^5(t 1 )q Lm , (5.9) 

where D tl = tijj-- Also we may assume that the summation in 15.91 goes from k = 1 3 . Or 
equivalently 

tv\MX(u, qtix\)X(v, Xi)q L( - ^ — tr|jvfX(ii, tiXi)X(v, Xi)q L ^ = 

Y^ c -^r^\ M {D k 5){ q n)q L ^\ 



Thus if we let F(h) = -r^- n tT\ M X(u 1 ,t 1 x 1 )X(v 1 ,x 1 )q L{0) , then 



chA/(<3) 

fc! 



F( 9 t 1 )-F(t 1 ) = ^^A fc 1 ^ 1 ) (5.10) 



If we multiply the equation above with (1 - ti) fel+1 we obtain (1 - ti) fe (F(gti) - F(ti)) = 0. By 
using Lemma 18 . II from Appendix A, it follows that the general solution of 15. 101 in the space 

q h C[[ti,ti 1 ,q\], 

is of the form 

F part (t 1 ) + f(q), 

where F par t(ti) is some particular solution and /(g) € C[[g]]. But F(ti) does not involve terms 
that contain only powers of q hence it follows that 



F(tl) = ££^ ( _ 1)fep f 1 



fe! 

Note that formally we have more than one solution of I5.1LW inside q h C[{ti, 1 , q, q^ 1 ]] (cf. 
Appendix A) . 

The 2-point functions are more interesting. 



tr\ M X(u 1 ,t 1 x 1 )X(v 1 , x^X (u 2 ,t 2 x 2 )X (v 2 , x 2 )q L< - 0) = 
tr\ M X(u!, ^—^-)X(v 1 ,x 1 )X(u 2 ,t 2 x 2 )X(v 2 ,x 2 )q L<C>) + 

q 

tT\ M [X(u 1 ,t 1 x 1 ),X(v 1 ,x 1 )X(u 2 , t 2 x 2 )X(v 2 ,x 2 )]q L< - 0) = 



i Note that for superalgebras this will not be the case 
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tr\MX(m, ^—^-)X(vi,Xl)X(U2,t 2 X 2 )X(v 2 ,X 2 )q L ^ + 

tr | M ^ CvigiM D%5(t 1 )X(u2,t 2 x 2 )X(v2,X2)q m + 
fcl>l 1- 



T ^r^tr\ M D^S (tm) X(v 1 ,x 1 )X(u2,X2t2)q LW = 
tr|MX(m, -!-^-)X(wi,a;i)X('U2,t2X2)X(u2,a;2)(7 i(0) + 

V ^^trj M A fc 1 1 5(ii)^(«2,t2^2)X(«2,x2) 9 L(0) + (5.11) 
fcl>l 

fc2>l o V / 



E^«^g(- 1 ) < ( i ?) J *" i '(^) 



X( ? ; 1 ,s 1 )^ 2 X(K2,tit2X 1 )g L(0) , 



fc 3 >i 

where fcj 6 N, i = 1, 2, 3. Since, 

trlwX^.xOD^-'XCiia.titaan)^ = trjw^-^^a, t -^-)X{v 1 , xi)q L(0) 

q 

we can extract zeroth term in 15. lit . For simplicity suppose that 114 = Vj = u, i = 1, 2 and let 
ki — k2 = ks = k. After taking the coefficient of x° in 15. 1H we obtain: 



If we act by q Dtl on 115. 121 1 we obtain 



F(qtx,t 2 ) - F(tx,t 2 ) = E ^ ( (9*0 F(t 2 ) + 

fe>i ' \ 

(-l) k D k 2 F (jjL^j + (-l) k D* 2 F( tl t 2 )). (5.13) 

By multiplying 15.131 with (1 — t\) K , where K is big enough, we obtain locality formula: 
{l-t x ) K {F(qtx,t 2 )-F(tx,t 2 )) = 

(l-«i) K E^((-l) k A^(^) +(-l) k D* 2 F(t 1 t 2 )). (5.14) 
The following proposition gives us an explicit formula for F{t\, t2). 
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Proposition 5.1 Suppose that 



H k {qtx,t 2 ) - H k (h,t 2 ) = (D 8) {qti)F(t 2 ) + 

t2_ 

qti 



DiF[^)+D^ 2 F(t 1 t 2 



(5.15) 



where F(t) = P k+1 ^ij . Then the 
is of the form 
(a) 



solution of \5. that satisfies Hk(ti,t 2 ) = Hk(t 2 ,ti) 



///.(/-i-/ >! =//('/)- l\ I i ( £,«J Pk+i 



(5.16) 



&) + If 



00 „2fc n 



(1 



+ £ 



n^g w ( (fafa) n +(tita)~ n ) 

(l-g~)2 



where fk{q) G C[[g]]. /n particular if /(g) = and |g| < 1 £/ien Hk(ti,t 2 ) converges 
uniformly on compact subsets inside a domain: 



\tit 2 q\ < 1, 



Ma 



< 1, 



< 1, 



t 2 



< 1, ti > 1, t 2 > 1. 



In particular, it has a meromorphic extension to (C x ) 2 
(b) IfU = e 2wiy -, i = l,2 andq = e 2lTiT , r£l, i/ien 

1 



Hi(t + yi,y 2 ) = H 1 (y u y 2 ) + 



and 



(2ni) 2 

H k {r + yx,y 2 ) = H k (y 1 ,y 2 



(pa(j/i - y 2 ,r) + p 2 (yi +y 2 ,r) + G 2 (r)) (5.17) 



(— 

{2m) k + 1 \ \dyi 



p 2 (yi - 2/2, t) + 



p 2 (2/i + 2/2,r 



(5.18) 



/or k > 2, where 

yi+y 2 i1®lT, 2/1 - y 2 £ Z © Zr. 

Proof: Proof of (a) requires solving 15.12L or I5.13H . inside q ,h C[[tf 1 , \ ?]]■ The Lemma 
gives us description of all solutions. Hence we have to find a particular one. This can be done in 
the following way (cf. [Sj, |2H])- We expand the left hand side of H5.12H in powers of t%. Then we 
compare coefficients by powers t™ on both sides, so we can determine coefficient of t\ in F(ti,t 2 ) 
(since it is multiplied by (1 — q~ n )). The details are not very illuminated so we omit them here. 
This gives us a way to find a particular solution of the g-difference equation 15.1211 . 

It is straightforward to see the uniform (and absolute) convergence of 15.16H when f(q) — 0. If 
is convergent. The locality formula 15.141 gives us a meromorphic continuation 
of lEHfl on (C x ) 2 such that for every n £ 1 

tit 2 / q n ,t\/t 2 / q n . 



Now if we apply the Proposition l5 . 1 1 we see that F(ti,t 2 ) that satisfies 15. 121 can be expressed 
as a linear combination in terms of Hk(ti,t 2 ), with fk(q) = 0. 
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Remark 5.2 Note that condition 12. 311 is crucial for our considerations. Otherwise one cannot 
obtain a recursion procedure, i.e. a way to express 4-point function by using 2-point functions. 
Zhu's recursion (cf. |29|~) would give us, in general, a way to express general 4-point function in 
terms of 3-point functions. 

Explicit calculation of the zeroth term of 15. Ill was done in [3j for V = M(l), in the case 
when Ui — Vi — h{— 1)1, i = 1, n, i.e., 

T(t 1 ,...,t n ,q)=7 1 (q)tr\ MW o(Y[h,t 1 ]h)---o(Y[h,t n }h)q Lm . (5.19) 

The approach in [3] is different since it deals with generalized characters and therefore explicit 
formulas are easier to obtain. In our approach we work with the formal delta function which 
carries an extra information not present in the generalized character formulation of 15.191 . 



6 Relation with vector bundles over elliptic curves 

Let Sir (t) be as before and E = C x /q z be an elliptic curve. Suppose that £(-q 1/2 t) G H {E,O* E ) 
is a holomorphic line bundle over E. Then £(— q 1,>2 t) has no global sections (the degree is -1). 
However F(t) = is a multi-valued section of £(~q 1/ ' 2 t). 

Suppose n = 2. Let t\ be the local coordinate on E and t 2 fixed. If we set 

Fi(ti) = G(*i,t 2 ) 

and 

F 2 (t 1 ) = G(t 1 t 2 ), 

where G(ti,ti) is Bloch-Okounkov 2-point function associated to a pair of fermions, then 

F^qh) = -g 1/2 tit 2 Fi(ti) + q 1/2 t 1 t 2 F 2 (t 2 ) 



F 2 (qt 1 ) = -g 1/2 ti£ 2 F 2 (£i) 



(6.1) 



Therefore a pair (Fi(ti), F 2 (t 2 )) can be though as a multi-valued section of the vector bundle V 
over E (of the rank 2) associated with the conjugacy class (which is enough to specify according 
to the classification theory) of the matrix 



- q 1/2 nt 2 





q 1/2 t lt2 
-q 1/2 t lt2 



The latter corresponds to a particular (non-split) extension of C{—q ' tit 2 ) by itself. 



This can be generalized for an arbitrary rank. Fix t 2 , 

l<il<...<ifc <n 



, t n . For 1 < k < n let 

i ? * * * ) ^ik J * " ')* 



(6.2) 



Notice that the right-hand side of 16.21 is invariant with respect to the symmetric group SVi-i 



(permuting t 2 ,...,t n ). Then (f\(ti) 
associated to 



, F„(ti)) is a multi-valued section of a rank n vector bundle 



1/2 



/ -i 





V - 



(6.3) 
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7 Conclusion and future work 



(a) Notice that one can easily generalize all results in our paper by considering more general 
iterated correlation functions of the form 

{u' n+1 , X(Y[ui,yi — Wi]vi, e wi x\) ■ ■ ■ X(Y[u„,y n — w n ]v n , e Wn x n )u n +i), 

{u' n+1 , X(ui,tixi)X(vi, sixi) ■ ■ ■ X(u n ,t„x H )X(v n , s n x n )u n+ i) , 
and the corresponding q-traces. 

(b) Our considerations in this paper were mostly at the level of vertex operator (super) algebras 
and its modules. If we closely examine the n-point function F(ti, ...,t„) associated to a 
pair of free charged fermions, we see that we are actually dealing with certain intertwining 
operators. More precisely, the space T§ — M(l) is a vertex operator subalgebra inside the 
vertex operator superalgebra T. Therefore operators X{ip,x) and X(t[>*,x) do not act on 
M(l) (even though the corresponding 2n-point trace can be computed for M(l) by using 
the boson-fermion correspondence) . Thus we are dealing with an non-trivial (even though 
simple) intertwining operators of the type 

where T±\ are M-modules. The fusion algebra for T, viewed as an intertwining operator 
algebra, is isomorphic to C[Z]. Therefore in the future we shall study iterated 2n-point 
correlation functions and corresponding q-traces for the abelian intertwining operator al- 
gebras. 

(c) It is possible to study, over the same lines, the n-point functions associated to certain 
N = 2 vertex operator superalgebras. 

(d) In I25| . motivated by we consider a large algebra of certain pseudodifferential opera- 
tors 2?* . In the Part III |2U we study n-point functions twisted with various Dirichlet 
characters. 

8 Appendix 

8.1 Appendix A 

In this appendix we prove some elementary results necessary for solving formal q-difference 
equations for certain 1- and 2-point functions. 

We treat the following two types of q-difference equations: 

F(qt)-q a F(t) =Bi(i) (8.1) 

and 

F(qti,ti) + q 1/2 tit 2 F(ti,t2) = B 2 {t u t 2 ), (8.2) 

whereB(t) G Cfc* 1 , t* 1 ]], B 2 (ti, t 2 ) G C^ 1 , tf 1 , tf 1 ]]. We discuss solutions both in C[[q ±x , tf 1 , t 
and C[[q, t^ 1 , tf 1 ]]. This is important since 

F(qt) -F(t) = D k t 8{qt) 
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has two distinguished solutions in C[[q ±:t , t^ 1 ]]. In general 18. Ill may not have any solution. For 
example 

F(qti) - F(t) = S , 
does not have solution in C[[q, tf 1 , tj 1 ]]. 

Remark 8.1 In most of the literature one does not consider formal solutions of q-difference 
equations but rather solutions in some analytic spaces (like holomorphic functions in the disk, 
punctured disk, etc.). We work often with series with no convergence (like delta functions) thus 
our approach has formal-distribution theoretical flavor. 

Lemma 8.1 Let ti£l Then all solutions of 

F(qt) - q n F(t) = 0, 

in CHq,^ 1 ]] are of the form t n f(q) for some f(q) £ C[[q]]. 

Lemma 8.2 All solutions of 

F(qt)+q 1/2 tF(t) = 0, 

in^'C^ 1 ,^ 1 ]] are of the form 

/(^(-lyV^r+l, (8.3) 

nGZ 

for some /(g) £ C[[g ±:L ]]. In particular the only solution inside t 1 / 2 C[[g, i ±:L ]] is the trivial one. 
Proof: Since F(t) = J2„ez fn+i/2(q)t n+ ^ , it follows that 

<Tfn+l/2(q) = -/n-l/2(<?)- 

Hence, /172(g) = -f-i/z(q) and fi/n(q) uniquely determine / n+1/2 (g). 
Lemma 8.3 All solutions of 

F(qti,t 2 )-F(t u t 2 )=0, 
inside C[[g ±:L , if 1 , t^ 1 ]] that satisfy F(ti,t2) = F(t 2 ,ti) are contained in C[[q ±:L ]]. 

Lemma 8.4 All solutions of 

F(qt 1 ,t 2 ) + q 1/2 t 1 t 2 F(t 1 ,t 2 ) = 0, (8.4) 
intl^tl^Cftq* 1 ,^ 1 ,^ 1 ]] that satisfy 

F(t 1 ,t 2 ) = F(t 2 ,ti) (8.5) 

are of the form 

/(<?)£ (-i)V^(ti*2) n+ i 

In particular there is no nontrivial solution inside C[[g, if 1 , i 2 *]]■ 
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Proof: Let F(ti,t 2 ) = X^nez+i ^1/1(^2, By comparing coefficients of t™ in 18,41 we get 

q n fn(t 2 ,q) +q 1/2 t 2 fn-l(t2,q) = 0, 
On the other hand 18.511 gives us 

fn{qt 2 ,q) +q 1/2 t 2 /n-i(i2,<?). 

Combined 

f„(qt 2 ,q)-q n f„(t2,q) = 0, (8.6) 

for every n £ Z + i. Because of Lemma l8.ll all solutions of 118.611 . inside i?/ 2 C[[q, t^ 1 ]], are of 
the form tj/(g), for some /(g) G Cfe* 1 ]]. Hence F(f u t a ) = £„ ez+ i (*i*a)"/n(«). Now apply 
Lemma l8.2l I 



8.2 Appendix B 

Here we give a different proof of the main technical result used by Zhu — the recursion formula 
(cf. Proposition 4.3.4 of [HO) 1 ) and its consequences. 

Instead of working with functions and their analytic properties 15.41 we prefer the use 

formal variable "all the way" and then, at the end, recognize certain formal series as Laurent 
expansions (around y = 0) of Weierstrass' functions. In other words, unlike Zhu's proof, our 
proof is completely formal — and then at the end — one can "turn on" the complex variables 

First we (re)prove the following result from 29 (Proposition 4.3.2): 

Proposition 8.1 

tr\ M X(u 1 ,x 1 )X(u 2 ,x 2 )q L( - 0) = trj A fo(ui)o(it 2 )(7 i(0) 

+ V P m+ i( — ,q)X( Ul [m]u 2 ,x 2 )q Lm . (8.7) 
c — ' X\ 

rn>0 

Proof: From |25). 

/ e v x 2 \ 

[X(ui,X!),X(u 2 ,x 2 )] = ReSyS I— — J X(Y[u!,y]u 2 ,x 2 ). 

Now from the property of the trace tiM^ABC) = trj m(BC A) we obtain 

(1 - q-°'i )X{u 1 ,x 1 )X{u 2 ,x 2 )q Lm = 

E (—) X( Ul [m]u2,x 2 ). (8.8) 

Now, by formally "inverting" the operator 1 — q~ Dlcl inside C[[q]] we obtain the formula 18.71 . I 

Let us recall Weierstrass function pk(z,r) defined before. These functions have a Laurent 
expansion near z = (k = 1, 2, ...) : 



n>l V 



G 2n+2 (T)z' n+ '~ k . (8.9) 
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Notice that 118.911 can be considered formally as an element of 

C[[2,2 -1 ,g]]. 

Let us define related function: 



- , . ) P2{z,q) + G 2 (q) for for k = 2 , . 

^ Z ^~< p k {z,q),k>Z (8 - 10) 



Theorem 8.1 



tr\ M X(Y[u,y]v,x) = 

£ p m +i{y)^\MX{u[m]v,x)q LW +tT\ M o{u)o(v)q L(0) . (8.11) 

m>l 

/n particular, if M is a V -module and V satisfies C 2 ~condition (cf. [29) then \H. lH is convergent 
and it has a meromorphic continuation to the whole y plane. 

Here are some consequences (cf. Propositions 4.3.4 and 4.3.5 in |30p: 
Corollary 8.1 

iT\ M o{u[-l]v)q L{0) = tr\ M o(u)o(v)q L(0) + £ G 2 k(q)tr\ M o(u[2k - 1]«) 9 L(0) . (8.12) 

m>l 

Combining Corollary 18.11 and Proposition 18 . 1 1 we obtain 
Corollary 8.2 

tr\ M X(u,x 1 )X(v,X2)q L(0) = (8.13) 
= tr\ M o(u[-l]v)q Lm + hr\ M o(u[0]v)q Lm + 

tr| M Yl P ™+! o(u[m]v)q Lm -^2tT\ M G 2k (q)o(u[2k- l]v)q Lm . 

m>0 ' fc>l 

Corollary 8.3 Let q = e 2wiT . The expression JO) is a doubly periodic (elliptic) function with 
respect to transformations 

y i— > y + 27ri, 
y i— > y + 27UT. 

Proof of Theorem \8. 1\ : 

tr\ M X(Y[u,y]v,X2)q LW = 

= Y, tr\ M X(u(i)v,x 2 )e yd ^ u \e y - 1)-^V (0) 

iSZ 

_ Jl eSx 2 ,deg(u) + dog(D)-i-l e ydcg(u)^ e y _ »-l 

iez 



tr| M ((an -x 2 ) i y( ?t ,x 1 )y(i;, : r 2 ) - (-x a + sti)^ («, x 2 )Y(u, an)) 9 L(0) 

= Res,, V e f dc sM( e f - l)-*- 1 ^ 1 ( - 

V : 



deg(u) 

-1 -1 12 ' 
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tr| M (x 2 l (xi -x 2 yX(u, x 1 )X(v,x 2 ) - x 2 \-x 2 + xi) % X(v, x 2 )X(u, xi) j g L(< " 

/ \ dcg(u) 

= Res xi Y e ydcg(u \e y - l) _i_1 Xa 1 — 



ifcS 



a; 2 l (xi - x 2 y y tr| M -Pm+i ( — , g ) o(u[m]v)q 
-x^\-X2 + x-i)" 1 y tr\ M P m +i f ?J o(w[m]w)g 



1,(0) 

deg(u) 



a; 2 l (a;i — £2)' — a; 2 £2 + xi) 1 tr\Mo(a)o(b)q 



M°) 



(8.14) 



We introduce a substitution i = 21 . Then 

a?2 l (a;i - x 2 y - x^. % {-xi + x ± y \ tv\ M o{a)o{b)q L{0) 



-- Rest £ e ydes(u) (e y - ip- V deg(,l) [ (t _ 1)* _ (_i + t f J tr| M o(a)o(6)g i(0) 

i€Z \ / 

= Res t ^ e ! ' dcs(u) (e ! '-l)- l - 1 t- dos(l,) (t-l) l tr|Mo(a) O (fe)g i(0) 

i<-l 

-Res t Yl (-i + tye ydce(u \e y -l)- l -h- dcs(u hr\ M o(a)o(b)q Lm 

»<-i 

= Rest E r^trl Mo(a)o(i)/' 0) 



i>0 



(t-l)< 

( „) (e^-1)* deg( „) , , . , . i(0 ) 



-Rest £ e» dcg <"> ^ - ^ ^ dcg( " ) tr| M o(a) (fe)g 

i>0 ^ ' 

= Rest ^'r^- 1 i _ 1 ey/t tr| M o(a)o(b)g £(0) ) + 
+Res t ^- d ^(") e ^(")- 1 T -L_tr| MO (a) O (6)g i(0) ) 

= tr| M o(a)o(6)g L(0) . (8.15) 
Also for every m > 1 

Re S:ci ^ e Hdcs(u) (e a - lP" 1 ^ 



. > deg(u) 
-1 -1 / ^2 



i£Z 
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X 2 '{Xl — X2) l Pm+l f ~~) QJ ~ x 2 \ — X2 + Xlf P m +1 {^T"^ 

= ResJ J2 e ydcs(u) (e y - lP -1 ^ 1 t" deg ("' • 
{(i - iyP m +l{t-\q) - {-l + t) l P m+ l{qr\q)} 

Res t { e »dcg( M ) t -dcg( M) -i _ _L_ p g) 



_|_^-deg(u) e J/deg(u)-l 



+Res J 1 t d ^ Y ( 4 ^-7^— jtt*(*) ■ (8.16) 

I e» - 1 ^ (e» - l) 1 (m-1)! w [ 



where we used the fact that 



P m (t,q) - P m {qt,q) = 7 — S(t). 



If we combine all P m 's into a single generating function we obtain 

y {Pm(t,q) - P m (gt,g))x™- 1 = 5(e x t). 



Now 



ydcg(u) 



Res f e " - " t dcg(u) Y ^ T^S(e x t) = 

e y _ 1 Z_/ ( e y _ J.)* v ' 

Rest- e -^s(u) y ye i> d ( e x t) 

e» -1 f-^ (e» - l) 1 

i=0 V ' 



e (y-a;)dog(u) 



(8.17) 



In order to evaluate the last sum in 118. 1611 we have to extract the coefficient of x m 1 inside 18.1711 . 
Clearly 

Coeff„ 



e y—x _ I 

= Coeff_. m -ie x d y — 

e y - I 

_ (-1)" 1 " 1 / d \ m_1 ^ B k (deg(u))y k ~ 1 



(m-1)! \dy J fc! 



y Cfc ™^ . (8.18) 



Now, by using the formula 



m — l „7ij.—n 



i I _ tti — 1 j.n 771 — 1 „ra ( 



(m-1)! 
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we obtain 



Res t { e !/dog( u ) f -do g ( u )-i _l_ Pm {r i^ q) + 

+r de g (u) e „de g C«)- 1 _l_p ro(3f -l )g) | = 

= (-1)"- 1 ^ cB ( dcg(u)+ „ ) g d ^")+"(deg( M )+n) m - 1 

(m - 1)! 1 - a deg(u)+n 

V ; n>0 y 

(m-1)!^ l-g n (m-1)! ^ 1 - g™ v ; 



From JCTfl and JO^ it follows (m > 2): 



Res,, y e^^e" _ l)- 1 - 1 ^ 1 — 

| ^a;^" ! (a;i — i2)'Pm g^j — % {— X2 + P m ^^~>9^ tr|Mo(u[m — l]w)g 
_/ 1 , (-l)" 1 " 1 V g m+ „(deg( U ))y" dcg ^" 1 



L(0) 



y m (m-1)! ^ n\ (m-1)! 

y v ' n>0 V ' n=l 



+7m^l)! E « WB W + (-D" 1 ^"^- I \^oHm - !]„), 



(m-1)! I ^ 1-g" 1-g 



i(0) 



1 , (-1)'" y / g m+r (deg( M )) 5 m+r \ y r | 
y m (m — 1)! 4^ \ m + r m + r J r\ 

E E r in-^ + rin-^ tr|M0(«[m - !]«)« 



(m-l)l Vj=o^i r!(l-g") 

_ f 1 (~l) m 

I 7/ m (m-1)! ^ \ (m + r)r! 

2 (r + m-l\ n m+r - 1 q n \ r \ . m 

— > y >tr \mo(u m — 1 \v)q y ' 

{r + m-iy.^y m-1 J 1-g" y y J J m 

_ f 1 (-l) m V 2 A n"-y \ 2fc -r 

tr|Mo(w[m— l]v)g i ^°' ) 

- f^ + (-l)" E f^I^G^^-^trlMoHm-lHg^) 
f (p 2 (?/,g) + G 2 (g))tr|Mo(i t [l]«)g i(0) , for m = 2 
(_ p m (t/,g)tr| M o(u[m - l]?;)g I ' (0) , for for m > 2. 
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= p m {y, g)tr| M o(u[m - l]v)q L(0) . (8.20) 
We used a fact that 

G ' 2k[T) = ~ (2fc)T (2^T)I ^ T^T- 

Because of 

trlMO^IO]^)^^ 5 = 0, (8.21) 
formulas jP2) and ftTST) imply 

tr| M ^(Vh,y]«,a;)g L(0) = 

pm+i(j/,g)tr| M X(w[m]i;,2:)g i(0) + tr| M o(u)o(u)g L(0) . (8.22) 

m>l 

Proof of Corollarv \8.,'Jl 

If we let g = e 2,rlT , the ellipticity with respect to 

yi-*-y + 'iitiT 

follows directly from formulas 15. II and 15.21 . I 
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